Gas Reservoir Evaluation and Assessment Tool 
5 Method and Apparatus and Program Storage Device 

Background of the Invention 

1 0 The subject matter of the present invention relates to computer modeling of a gas 

reservoir, and, more particularly, to a method and apparatus and program storage device 
adapted for generating a computer model which will predict the pressure and the 
production behavior of a gas reservoir. 

1 5 Reservoir simulation is an essential tool for the management of oil and gas reservoirs. 

Prediction of the pressure in and the production of a gas reservoir under various operating 
conditions allows, among other benefits, proper investment decisions to be made. In 
order to make such a prediction, one must construct a reservoir model. The reservoir 
model is essentially a mathematical model that is implemented through a computer 

20 program. History matching observed behavior of the reservoir must validate the 
parameters of the model. Ideally, finite difference simulators are used to construct 
reservoir models. This permits detailed characterization including heterogeneity, 
multiphase effects like water coning and fingering. However, in order to make full use of 
such a tool, a large amount of reliable data is required. Also, a full study, including the 

25 history matching step, may take months to carry out., Therefore, there is a demand for an 
alternative tool that honors the physics of fluid flow and, at the same time, generates a 
solution which is many orders quicker and faster than the aforementioned finite difference 
simulator. 
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Summary of the Invention 



One aspect of the present invention involves a method of generating a prediction of 
values in a reservoir, comprising the steps of: (a) receiving input data characterizing the 
reservoir; (b) producing a computer model in response to the input data representing the 
reservoir, the producing step (b) of producing the computer model including the steps of, 
(bl) calculating the values in one dimension associated with a single layer in the 
reservoir, each of the values existing at a single point in space in the reservoir and at a 
single point in time in the reservoir, (b2) calculating the values in the one dimension 
associated with multiple layers in the reservoir, each of the values in each of the multiple 
layers existing at a single point in space in the reservoir and at a single point in time in the 
reservoir, (b3) calculating the values in three dimensions associated with the multiple 
layers in the reservoir, each of the values in each of the multiple layers in the three 
dimensions existing at a single point in space in the reservoir and at a single point in time 
in the reservoir, (b4) calculating the values in the three dimensions as a function of time, 
the values being associated with the multiple layers in the reservoir, each of the values in 
each of the multiple layers in the three dimensions existing at a single point in space in 
the reservoir, the each of the values in the each of the multiple layers in the three 
dimensions existing at any future point in time in the reservoir, the computer model being 
produced in response to the calculating step (b4); verifying the computer model; and 
using the computer model, generating the prediction of the values in the reservoir in 
response to the verifying step. 

Another aspect of the present invention involves a program storage device readable by a 
machine tangibly embodying a set of instructions executable by the machine to perform 
method steps for generating a prediction of values in a reservoir, the method steps 
comprising: (a) receiving input data characterizing the reservoir; (b) producing a 
computer model in response to the input data representing the reservoir, the producing 
step (b) of producing the computer model including the steps of, (bl) calculating the 



values in one dimension associated with a single layer in the reservoir, each of the values 
existing at a single point in space in the reservoir and at a single point in time in the 
reservoir, (b2) calculating the values in the one dimension associated with multiple layers 
in the reservoir, each of the values in each of the multiple layers existing at a single point 
in space in the reservoir and at a single point in time in the reservoir, (b3) calculating the 
values in three dimensions associated with the multiple layers in the reservoir, each of the 
values in each of the multiple layers in the three dimensions existing at a single point in 
space in the reservoir and at a single point in time in the reservoir, (b4) calculating the 
values in the three dimensions as a function of time, the values being associated with the 
multiple layers in the reservoir, each of the values in each of the multiple layers in the 
three dimensions existing at a single point in space in the reservoir, the each of the values 
in the each of the multiple layers in the three dimensions existing at any future point in 
time in the reservoir, the computer model being produced in response to the calculating 
step (b4); verifying the computer model; and using the computer model, generating the 
prediction of the values in the reservoir in response to the verifying step. 

Another aspect of the present invention involves a system adapted for generating a 
prediction of values in a reservoir, comprising: first apparatus adapted for receiving input 
data characterizing the reservoir; second apparatus adapted for producing a computer 
model in response to the input data representing the reservoir, the second apparatus 
adapted for producing the computer model including, third apparatus adapted for 
calculating the values in one dimension associated with a single layer in the reservoir, 
each of the values existing at a single point in space in the reservoir and at a single point 
in time in the reservoir, fourth apparatus adapted for calculating the values in the one 
dimension associated with multiple layers in the reservoir, each of the values in each of 
the multiple layers existing at a single point in space in the reservoir and at a single point 
in time in the reservoir, fifth apparatus adapted for calculating the values in three 
dimensions associated with the multiple layers in the reservoir, each of the values in each 
of the multiple layers in the three dimensions existing at a single point in space in the 
reservoir and at a single point in time in the reservoir, sixth apparatus adapted for 



calculating the values in the three dimensions as a function of time, the values being 
associated with the multiple layers in the reservoir, each of the values in each of the 
multiple layers in the three dimensions existing at a single point in space in the reservoir, 
the each of the values in the each of the multiple layers in the three dimensions existing at 
5 any future point in time in the reservoir, the computer model being produced in response 
to the calculating performed by the sixth apparatus; seventh apparatus adapted for 
verifying the computer model thereby generating a verified computer model; and eighth 
apparatus, responsive to the verified computer model, adapted for generating the 
prediction of the values in the reservoir in response to the verifying performed by the 
10 seventh apparatus. 

Another aspect of the present invention involves a method of producing a computer 
model in response to input data representing a reservoir, comprising the steps of: (a) 
calculating values in one dimension associated with a single layer in the reservoir, each of 

1 5 the values existing at a single point in space in the reservoir and at a single point in time 
in the reservoir, (b) calculating the values in the one dimension associated with multiple 
layers in the reservoir, each of the values in each of the multiple layers existing at a single 
point in space in the reservoir and at a single point in time in the reservoir, (c) calculating 
the values in three dimensions associated with the multiple layers in the reservoir, each of 

20 the values in each of the multiple layers in the three dimensions existing at a single point 
in space in the reservoir and at a single point in time in the reservoir, and (d) calculating 
the values in the three dimensions as a function of time, the values being associated with 
the multiple layers in the reservoir, each of the values in each of the multiple layers in the 
three dimensions existing at a single point in space in the reservoir, the each of the values 

15 in the each of the multiple layers in the three dimensions existing at any future point in 
time in the reservoir, the computer model being produced in response to the calculating 
step (d). 
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Another aspect of the present invention involves a program storage device readable by a machine 
tangibly embodying a program of instructions executable by the machine to perform method steps for 
producing a computer model in response to input data representing a reservoir, said method steps 
5 comprising: (a) calculating values in one dimension associated with a single layer in said reservoir, 
each of said values existing at a single point in space in said reservoir and at a single point in time in 
said reservoir, (b) calculating said values in said one dimension associated with multiple layers in said 
reservoir, each of said values in each of said multiple layers existing at a single point in space in said 
reservoir and at a single point in time in said reservoir, (c) calculating said values in three dimensions 

10 associated with said multiple layers in said reservoir, each of said values in each of said multiple 

layers in said three dimensions existing at a single point in space in said reservoir and at a single point 
in time in said reservoir, and (d) calculating said values in said three dimensions as a function of time, 
said values being associated with said multiple layers in said reservoir, each of said values in each of 
said multiple layers in said three dimensions existing at a single point in space in said reservoir, said 

1 5 each of said values in said each of said multiple layers in said three dimensions existing at any future 
point in time in said reservoir, said computer model being produced in response to the calculating 
step (d). 

Another aspect of the present invention involves a system adapted for producing a computer model in 
20 response to input data representing a reservoir, comprising: first apparatus adapted for calculating 

values in one dimension associated with a single layer in said reservoir, each of said values existing at 
a single point in space in said reservoir and at a single point in time in said reservoir, second 
apparatus adapted for calculating said values in said one dimension associated with multiple layers in 
said reservoir, each of said values in each of said multiple layers existing at a single point in space in 
25 said reservoir and at a single point in time in said reservoir, third apparatus adapted for calculating 
said values in three dimensions associated with said multiple layers in said reservoir, each of said 
values in each of said multiple layers in said three dimensions existing at a single point in space in 
said reservoir and at a single point in time in said reservoir, and fourth apparatus adapted for 
calculating said values in said three dimensions as a function of time, said values being associated 
30 with said multiple layers in said reservoir, each of said values in each of said multiple layers in said 
three dimensions existing at a single point in space in said reservoir, said each of said values in said 
each of said multiple layers in said three dimensions existing at any future point in time in said 
reservoir, said computer model being produced when said fourth apparatus calculates said values in 
said three dimensions as a function of time. 
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Brief Description of the Drawings 



A full understanding of the present invention will be obtained from the detailed description of the 
5 preferred embodiment presented hereinbelow, and the accompanying drawings, which are given 
by way of illustration only and are not intended to be limitative of the present invention, and 
wherein: 



figures 1 and 2 illustrate a method of computer modeling for predicting the pressure in a gas 
reservoir and the production of gas from the gas reservoir using a finite difference engine and 
analytical engine, respectively; 



figure 3 illustrates a workstation which is responsive to certain specific 'input data' and which 
stores a Gas Reservoir Evaluation and Assessment software of a Gas Reservoir Evaluation and 
1 5 Assessment Tool (GREAT) in accordance with the present invention; 

figures 4 and 5 illustrate a block diagram of the Gas Reservoir Evaluation and Assessment 
software of figure 3 which includes an 'Analytical Engine', and the 'input data' which is 
introduced to the Gas Reservoir Evaluation and Assessment software of figure 3; 

20 

figures 6 and 7 illustrate a construction of the 'verification' step and the 'Analytical Engine' step, 
respectively, of the Gas Reservoir Evaluation and Assessment software of figure 4 of the present 
invention; 

15 figures 8-1 1 illustrate cross sectional views of an Earth formation designed to assist in an 
understanding of figure 7 of the drawings which represents an illustration of the Analytical 
Engine; and 



figure 12 illustrates an example of the Output Record which is generated by the Recorder or 
Display Device associated with the Gas Reservoir Evaluation and Assessment Tool (GREAT) of 
figure 3. 
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Description of the Invention 



Prediction of the pressure-production behavior of a hydrocarbon reservoir is essential for its 
5 effective management. Project planning and screening functions depend on time 

availability of such information. There is a need for a fast solution that involves history 
matching and subsequent prediction. A Gas Reservoir Evaluation and Assessment Tool 
(GREAT) disclosed in this specification is based around a newly formulated set of 
equations applicable to multiple wells in a single-phase system. The GREAT tool provides 

1 0 a complete workflow for gas reservoir evaluation comprised of data entry and 

manipulation, model initialization, test interpretation, history matching and prediction. The 
GREAT tool includes the Analytical Engine 20 which further includes a newly derived 
solution of diffusivity equations for multiple wells, horizontal or vertical, in a single phase 
layered system under a variety of boundary conditions. The solution of these equations 

1 5 model both transient and steady state flow regimes and is applicable to both testing and 
long term performance prediction. The equations applicable to laminar flow of fluids in a 
porous medium were the results of Darcy's experimental study of the flow characteristics of 
sand filters. This combined with the equation of continuity and an equation of state for 
slightly compressible fluid yields the diffusivity equation, which is the equation for 

20 pressure diffusion in a porous medium. Solution of the diffusivity equation under different 
boundary conditions forms the basis for a prediction of the bottom hole pressure response 
of a producing well. These analytical solutions are generally applicable for a single well 
and are used widely in the area of well testing. The efficiency of analytical models is 
generally judged by accuracy and speed. The novel set of solutions used in the GREAT 

25 tool is applicable to multiple wells, which can be vertical as well as horizontal. These wells 
can be operating as producers or injectors thus being of additional significance to gas well 
storage. The solution of the diffusivity equation set forth in this specification has been 
derived by application of successive integral transforms. The application of these new 
solutions is characterized by stability and speed. 
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Accordingly, in this specification, a Gas Reservoir Evaluation and Assessment Tool 
(GREAT) in accordance with the present invention utilizes an Analytical Engine (instead 
of a Finite Difference Engine) to produce predictions of pressure values and other 
production data at 'any point in space' and at 'any point in time' in a reservoir. A 
computer system, such as a workstation, stores a Gas Reservoir Evaluation and 
Assessment software which includes the Analytical Engine and responds to input data 
(which includes a reservoir description and fluid properties) by generating an output 
record which represents a prediction of the pressure values and other data at 'any point in 
space' and at 'any point in time' in a reservoir. The Analytical Engine will first calculate 
a pressure value in ID for a single layer of a reservoir at a 'single point in space' and a 
'single point in time'; it will then calculate a pressure value in ID for multiple layers in 
the reservoir at the 'single point in space' and the 'single point in time'; it will then 
calculate a pressure value in 2D for the multiple layers at the 'single point in space' and 
the 'single point in time'; it will then calculate a pressure value in 3D for the multiple 
layers at the 'single point in space' and the 'single point in time'; and it will then 
calculate a pressure value in 3D for multiple layers not only at a 'single point in space' 
but also at 'any future point in time'. 

Referring to figure 1, a first method of computer modeling is illustrated adapted for 
predicting the pressure in a gas reservoir and/or the production of gas from the gas 
reservoir, this first method using a 'finite difference engine'. In figure 1, the first step of 
the method is Data Entry step 1 0. This Data Entry step 1 0 involves the entry of input 
data representing a Gas Reservoir, such as Well Radius, skin which represents a zone of 
damage around a wellbore, perforation length, fluid properties, a reservoir description, 
and rate data. This input data will be discussed in greater detail later in this specification. 
A 'Finite Difference Engine' 12 receives the input data provided via the Data Entry step 
10. The 'Finite Difference Engine' 12 includes a 'Model Initialization' step 12a. 
Usually, the 'Model Initialization' step 12a includes a Pressure Calculator for calculating 
a pressure in the Gas Reservoir at a single point in space and at a single point in time 
(hereinafter called a 'pressure in ID at a point in time', where ID means 'one dimension). 
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The 'Model Initialization 5 step 12a will also be discussed later in this specification. The 
'Finite Difference Engine' 12 is basically a 'numerical simulation method' where the 
simulation or modeling is performed 'numerically'. The 'Finite Difference Engine' 12 is 
error prone since errors can accumulate during the numerical simulation process. In 
5 addition, the numerical simulation process associated with the 'Finite Difference Engine' 
12 is very time consuming since a large amount of time is consumed during the numerical 
simulation process. The 'Finite Difference Engine' 12 will produce a 'computer model' 
which must be verified during the Verification step 14. During the Verification step 14, 
'known data' (having 'known results') will be introduced into the 'computer model' that 

10 was received from the 'Finite Difference Engine' 12. Responsive thereto, the 'computer 
model' will generate 'results'. Those 'results' are compared with the aforementioned 
previously 'known results'. If the 'results' approximately equal the 'known results', the 
Verification step 14 reveals that the 'computer model' which was generated from the 
'Finite Difference Engine' 12 does, in fact, produce accurate results and accurate 

1 5 predictions. The Verification step 1 4 includes two types of verification: ( 1 ) Test 

Interpretation 14a, and (2) History Matching 14b (each verification step 14a and 14b will 
be described later in greater detail). Following the Verification step 14, since the 
'computer model' generated by the 'Finite Difference Engine 12 was properly verified 
(i.e., the 'computer model' does, in fact, produce accurate results or predictions), the next 

20 step includes the Prediction step 16. At this point, the 'computer model', generated by 
the 'Finite Difference Engine' 12 and verified during the Verification step 14, can predict 
the pressure in, and/or the production from, a Gas Reservoir at any point in space in the 
reservoir and at any future point in time in the reservoir. As a result, during the 
Prediction step 16, a time variable 't', representing a 'future point in time t', can be 

25 entered during the 'Data Entry' step 10, the 'enter' key on a computer system or 

workstation is depressed, and, as a result, the 'pressure' in the Gas Reservoir and/or other 
production data associated with the Gas Reservoir, for that particular future point in time 
't', will be displayed or otherwise provided to a user. However, as previously mentioned, 
during the 'Finite Difference Engine' step 12, the 'Finite Difference Engine' 12 

30 processing time can be very 'time consuming' and, in addition, the results produced by 
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the 'Finite Difference Engine' 12 can include a large quantity of errors, since a number of 
individual errors can accumulate during the 'Finite Difference Engine' stepl2. 

Referring to figure 2, a second method of computer modeling is illustrated adapted for 
5 predicting the pressure in a gas reservoir and/or the production of gas from the gas 
reservoir, this second method using an 'analytical engine 5 . In figure 2, this second 
method is identical to the first method of figure 1, except that the second method of figure 
2 uses an 'Analytical Engine 5 which replaces the 'Finite Difference Engine' of figure 1. 
In figure 2, the Data Entry step 10 is followed by the 'Analytical Engine' step 20. The 

10 'Analytical Engine' step 20 includes the 'Model Initialization' step 12a as in figure 1. 
The 'Analytical Engine' step 20 is followed by the same Verification step 14 as in figure 
1 and the Verification step 14 is followed by a 'Prediction' step 18. The 'Prediction' step 
18 is not the same as the 'Prediction' step 16 of figure 1 because the 'prediction' or 
'results' 18 provided by the 'Prediction' step 18 are much more accurate and, in addition, 

1 5 that 'prediction' 1 8 is provided instantaneously as noted below. As mentioned earlier, 
the 'Analytical Engine' 20 is orders of magnitude faster than the 'Finite Difference 
Engine' 12 of figure 1 because the 'Analytical Engine' 20 provides the 'results' or 
'prediction' 18 nearly instantaneously following the Data Entry step 10 of figure 2. In 
fact, when the Data Entry step 10 of figure 2 is complete, and when the 'enter' key on the 

20 workstation or computer keyboard is depressed, the 'prediction' from the 'prediction' 
step 18 is instantaneously available to a user via a Recorder or Display device of a 
workstation or computer system. In addition, the 'prediction' 18 produced by the 
'Analytical Engine 5 20 is much more accurate than the 'prediction' 16 produced by the 
'Finite Difference Engine' 12. A discussion of the 'Analytical Engine' 20 will be set 

25 forth below in the 'Description of the Invention' portion of the specification; however, in 
addition, a detailed specification document is provided in the 'Detailed Description of the 
Invention' portion of this specification, that detailed specification document providing a 
detailed construction of the 'Analytical Engine' 20 of figure 2. 

30 In figures 1 and 2, the Data Entry step 10 relates to entry of historical rate and pressure 

10 



data. Pressure data may include static or flowing bottom hole pressure measurements of 
each well over the producing history of the reservoir. In addition, it could also include 
pressure transient data. The GREAT tool of the present invention shown in figure 3 will 
be particularly suitable for handling data from permanent downhole gauges. Permanent 
5 downhole gauges provide a large quantity of data captured at very small time intervals 
over extended periods of time. Suitable filtering techniques will be provided to process 
this data for subsequent use. Processed well logs can be entered and displayed. 

The Model Initialization step 12a will allow for initiation of the gas reservoir model with 
10 a basic model shape, layer geometry, and static properties, such as porosity. A range of 
correlation suites will be provided to generate fluid properties. 

During the Test Interpretation step 14a, initial estimates of permeability and skin, 
D- factor, open flow potential, initial reservoir pressure, and reservoir volume can be 
1 5 made. Basic line fitting functionality on specialized transient pressure plots will be used 
for this purpose. Additionally, two basic material balance plots, i.e., Havlena - Odeh and 
p/Z, will be available to determine the reservoir volume from static pressure data. 

The History Matching step 14b is essentially for model validation. Having obtained the 
20 initial parameter estimates from previous steps, non-linear regression can now be 

performed on the observed pressure data for all the wells in the reservoir. This would 
provide both final tuning of the reservoir and well parameters in the context of the entire 
reservoir. 

25 The Prediction step 16 and 18 will use the tuned model to provide forecast of reservoir 
and well bottom hole pressures based on individual well targets. Individual well targets 
can be set from Daily Contracted Quantity (DCQ) and monthly profile factors. Individual 
well and reservoir pressure cutoffs along with swing factor will provide the necessary 
checks for determining the ability of a well to produce. This step will also provide an 

30 option to generate a numerical model of the reservoir for simulation engineers. 
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It should be noted that it is not essential to follow all the steps. If pressure transient data 
is not available, initial permeability values can be obtained from logs and reservoir 
volume can be obtained from geological information. One can move from the Model 
Initialization step to the History Matching step. Similarly, if no historical data is present, 
one may use the Prediction step as a design exercise. 

Referring to figure 3, a workstation 22 or other computer system 22, representing a Gas 
Reservoir Evaluation and Assessment Tool (GREAT), is illustrated. In figure 3, the 
workstation 22 includes a processor 22a operatively connected to a system bus, a 
Recorder or Display Device 22b operatively connected to the system bus, and a memory 
or other program storage device 22c operatively connected to the system bus. The 
program storage device 22c is adapted to receive a 'Gas Reservoir Evaluation and 
Assessment software' 22c 1 from a CD-Rom or other such program storage device, the 
'Gas Reservoir Evaluation and Assessment software' 22c 1 being loaded from the CD- 
Rom into the memory/program storage device 22c of figure 3 for storage therein. The 
system bus of the workstation 22 is adapted to receive certain specific 'input data' 24, the 
'input data' 24 being discussed below with reference to figure 5. When the processor 22a 
of the workstation 22 receives the 'input data' 24 and executes the 'Gas Reservoir 
Evaluation and Assessment software' 22c 1 stored in the memory/program storage device 
22c using such 'input data' 24, an 'Output Record' 22bl is recorded or displayed on the 
Recorder or Display device 22b of the workstation or other such computer system 22 of 
figure 3. That 'Output Record' 22bl is discussed below with reference to figure 12. The 
workstation 22 or other such computer system 22 may be a personal computer (PC), a 
workstation, or a mainframe. Examples of possible workstations include a Silicon 
Graphics Indigo 2 workstation or a Sun SPARC workstation or a Sun ULTRA 
workstation or a Sun BLADE workstation. The memory/program storage device 22c or 
the aforementioned CD-Rom is a computer readable medium or a program storage device 
which is readable by a machine, such as the processor 22a. The processor 22a may be, 
for example, a microprocessor, microcontroller, or a mainframe or workstation processor. 
The memory/program storage device 22c, which stores the Gas Reservoir Evaluation and 
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Assessment software 22c 1, may be, for example, a hard disk, ROM, CD-ROM, DRAM, 
or other RAM, flash memory, magnetic storage, optical storage, registers, or other 
volatile and/or non-volatile memory. 

Referring to figure 4, the Gas Reservoir Evaluation and Assessment software 22c 1 of 
figure 3 includes four major steps, which steps are also illustrated in figure 2: (1) a Data 
Entry step 10, (2) an Analytical Engine step 20, (3) a Verification step 14, and (4) a 
Prediction step 18. The 'Output Record' 22bl of figure 3 is generated following the 
Prediction step 18. Each step is discussed below, as follows. 

Referring to figures 3, 4 and 5, during the Data Entry step 10 of figure 4, the 'Input Data' 
24 of figure 3 is entered into the workstation 22 of figure 3 via a workstation keyboard. 
Recall that the Gas Reservoir Evaluation and Assessment software 22c 1 is adapted to 
evaluate a gas (or oil) reservoir, where the reservoir includes a plurality of perforated 
wellbores penetrating an Earth formation, each wellbore having a well radius and a skin 
on a wall of the wellbore for producing a fluid from the formation. The gas (or oil) 
reservoir has a 'reservoir description' which includes 'reservoir data' that characterizes 
the formation penetrated by the one or more wellbores in the gas (or oil) reservoir, the 
'reservoir data' including permeability, resistivity, porosity, and the presence or absence 
of boundaries in the wellbore which may impede the flow of the fluid from the formation, 
etc. As a result, in figure 5, during the Data Entry step 10 of figure 4, 'Input Data' 24 will 
be entered into the workstation 22 of figure 3, where that 'Input Data' 24 includes: (1) the 
'well radius' representing a radius of a wellbore in the reservoir, (2) the 'skin' 
representing a 'zone of damage' surrounding the wellbore where the damage results from 
perforations in the formation penetrated by wellbore and from drilling fluid absorbed by 
the formation which reduces the ability of oil or gas to flow from the formation, (3) 
'perforation length' representing the length of perforations in the formation penetrated by 
the wellbore, (4) 'fluid properties' representing properties of a fluid produced by the 
wellbore, such as density, and other fluid properties, such as the pressure and temperature 
at which gas breaks out of solution, (5) 'reservoir description' representing data that 
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characterizes the formation penetrated by wellbores in the reservoir, such as the 
permeability, resistivity, or porosity of the formation, and the presence or absence of 
boundaries in the formation penetrated by the wellbore which may impede the flow of gas 
or fluid out of the perforations in the formation surrounding the wellbore, and (6) 'rate 
data' representing the rate at which gas (or oil) is produced at the surface of the wellbore 
as a result of the production of gas (or oil) from the perforations in the formation 
penetrated by the wellbore. 

Referring to figure 4, during the Analytical Engine step 20, a 'computer model' 20a will 
be generated. That 'computer model' 20a is adapted to predict the pressure of a gas in the 
gas reservoir and the production characteristics of the gas from the gas reservoir at any 
point in space and at any point in time. The Analytical Engine step 20 will be discussed 
below in greater detail in the 'Description of the Invention' section of this specification 
with reference to figure 7 and in the 'Detailed Description of the Invention' section of this 
specification where a specification document is set forth describing in detail the 
characteristics of the Analytical Engine 20. The Analytical Engine 20 produces the 
'computer model' 20a almost instantaneously following the Data Entry step 10 (when the 
'enter' key is depressed on a keyboard of workstation 22) and the 'computer model' 20a 
will produce predictions 18 in figures 2 and 4 which are much more accurate than are the 
predictions 16 of figure 1. 

Referring to figures 4 and 6, during the Verification step 14, the 'computer model' 20a is 
verified to ensure that the predictions' 18 will be accurate. For example, a pressure gauge 
will be lowered downhole into a wellbore at a 'particular point in space' and at a 
'particular point in time' to record a first pressure in the wellbore. Then, the 'computer 
model' 20a is interrogated responsive to the 'particular point in space' and the 'particular 
point in time' to produce a second pressure. If the first pressure is approximately equal to 
the second pressure, the 'computer model' 20a has been 'verified' to produce accurate 
predictions 18. In figure 6, the Verification step 14 includes two types of verification: (1) 
Test Interpretation 14a, and (2) History Matching 14b. In Test Interpretation 14a, 
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compare a prediction of a pressure in a reservoir as a function of space and time against 
simpler known cases. In History Matching 14b, compare a prediction of pressure in a 
reservoir as a function of space and time against a known measurement as a function of 
that same space and time. The time scales and the methodologies used during the Test 
5 Interpretation step 14a and the History Matching step 14b are different, but they both 
have the same purpose; that is, to verify that the 'computer model' 20a will produce 
predictions 18 that are accurate, as described above. 

Referring to figures 3 and 4, during the Prediction step 18 of figure 4, now that the 
10 'computer model' 20a has been determined to produce accurate predictions 1 8, the 
'computer model' 20a can be interrogated to produce a 'real time' prediction 18. 
Recalling that the 'computer model' 20a represents the reservoir whose characteristics are 
set forth in the Input Data 24 of figure 5, two types of data are provided to the 'computer 
model' 20a: (1) a 'future time' representing a point in time in the reservoir when the 
15 characteristics of the reservoir are desired, and (2) a 'point in space' representing a 
specific location in the reservoir where the characteristics of the reservoir are desired. 
Responsive thereto, the 'computer model' 20a will generate a 'Prediction' 18 which 
includes the 'pressure' and 'other production characteristics' in the reservoir at that 'point 
in space' and at that 'future time 5 . This 'pressure' and 'other production characteristics' 
20 information is valuable to a customer because the customer will want to know how the 
reservoir will produce in five years, for example, or when the reservoir will be depleted, 
etc. When the Prediction step 18 is completed, an 'Output Record' 22bl is generated 
which reflects and/or records and/or displays that prediction 18, that 'Output Record' 
22bl being generated by the Recorder or Display device 22b of figure 3. 

25 

Referring to figures 7 through 1 1, a construction of the Analytical Engine 20 of figures 2 
and 4 is illustrated in figure 7. 

In figure 7, the Analytical Engine 20 includes a Model Initialization step 12a which 
30 responds to the Input Data 24 of figures 3 and 5 including the well radius, skin, 
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perforation length, fluid properties, reservoir description, and rate data. In figure 7, the 
Model Initialization step 12a is a 'pressure calculator' 12a which calculates a pressure in 
a reservoir as a function of space and time; that is, the 'pressure calculator' 12a will 
calculate a plurality of pressures in one-dimension (ID), such as along the x-axis in the 
reservoir, where each of the plurality of pressures along the one-dimensional x-axis in the 
reservoir represent a pressure at a single point in space and a single point in time in the 
reservoir. Since the 'pressure calculator' 12a of the Model Initialization step 12a 
calculates a plurality of pressures in ID along the x-axis in the reservoir (each pressure 
representing a pressure at a single point in space and a single point in time in the 
reservoir), the 'pressure calculator' 12a actually calculates the plurality of pressures by 
performing a summation in the form of an 'Integral from 0 to x' along the reservoir, as 
follows: £ 

In figure 8, for example, a cross sectional view of an Earth formation is illustrated. 
Figure 8 represents a single layer of Earth formation. A plurality of pressure values 28 
have been calculated, during the Model Initialization step 12a of figure 7, along the x-axis 
in the single layer of Earth formation, each of the plurality of pressure values 28 
representing a pressure which exists at a 'single point in space' and a 'single point in 
time' in the reservoir. 

In figure 7, the Analytical Engine 20 also includes a Multi-Layer Loop step 26. The 
Model Initialization step 12a assumed that a 'single layer of formation' existed in the 
Earth formation of the reservoir when the plurality of pressure values 28 were calculated 
in a 'single point in space' within the 'single layer of formation' and during a 'single 
point in time' . However, we all know that the Earth formation consists of a multitude of 
layers (and not merely of a single layer as illustrated in figure 8). Therefore, the Multi- 
Layer Loop step 26 of figure 7 will respond to the 'Reservoir Description' portion of the 
Input Data 24 by converting the 'pressure in ID for one layer of the reservoir' (that was 
generated by the 'pressure calculator' 12a of the Model Initialization step 12a) into a 
'multi-layer formation' which is based on the reservoir description, the 'multi-layer 
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formation' representing a 'plurality of Earth formation layers'. In each of the 'plurality of 
Earth formation layers', a plurality of pressure values are calculated in ID along the x- 
axis of each layer, each pressure value being located at a 'single point in space' and 
existing at a 'single point in time'. 

In figure 9, for example, whereas in figure 8 a single layer of Earth formation existed, in 
figure 9, the Multi-Layer Loop 26 of figure 7 responds to the 'reservoir description' by 
converting the 'pressure in ID for one layer of the reservoir' that was generated by the 
pressure calculator 12a during the Model Initialization step 12a into a 'plurality of Earth 
formation layers' 36a, 36b, and 36c in the Earth formation. In figure 9, in each of the 
'plurality of Earth formation layers' 36a-36c generated by the Multi-Layer Loop 26 of 
figure 7, a 'plurality of pressure values' have been calculated along the x-axis of the 
formation at a 'single point in space' and at a 'single point in time'. For example, in 
figure 9, a 'plurality of pressure values' 30 have each been calculated in one-dimension 
(ID) along the x-axis of the formation layer at a 'single point in space' and a 'single point 
in time' for Earth formation layer 36a. In figure 9, a 'plurality of pressure values' 32 
have each been calculated in ID along the x-axis at a 'single point in space' and a 'single 
point in time' for Earth formation layer 36b. In figure 9, a 'plurality of pressure values' 
34 have each been calculated in ID along the x-axis at a 'single point in space' and a 
'single point in time' for Earth formation layer 36c. 

In figure 7, in connection with the Pressure Source Integrand 38, recall from figure 9 that 
a 'plurality of pressure values' 30, 32, and 34 have each been calculated in ID at a 'single 
point in space' and a 'single point in time' for the plurality of formation layers 36a, 36b, 
and 36c. In figures 7 and 10, the Pressure Source Integrand 38 of figure 7 will add an 
additional dimension to each formation layer 36a, 36b, 36c by converting the 'one 
dimensional (ID)' calculation of pressure values 30-34 for each formation layer 36a, 36b, 
and 36c (at the 'single point in space' and the 'single point in time') as shown in figure 9 
into a 'two dimensional (2D)' calculation of pressure values for each formation layer 36a, 
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36b, and 36c (at the 'single point in space' and the 'single point in time') as shown in 
figure 10. 

In figure 10, for example, the Pressure Source Integrand 38 of figure 7 will calculate an 
'additional plurality of pressure values' 40, 42, 44 along the y-axis of each formation 
layer 36a, 36b, 36c, each pressure value of the 'additional plurality of pressure values' 40, 
42, 44 being calculated at a 'single point in space' and during a 'single point in time'. 
Therefore, in figure 10, the 'plurality of pressure values' 30, 32, 34 have each been 
calculated in connection with the x-axis in each formation layer 36a, 36b, 36c at a 'single 
point in space' and during a 'single point in time'; and, in addition, the 'additional 
plurality of pressure values' 40, 42, 44 have each been calculated in connection with the 
y-axis in each formation layer 36a, 36b, and 36c at a 'single point in space' and during a 
'single point in time'. The Pressure Source Integrand 38 of figure 7 actually calculates 
the 'additional plurality of pressure values' 40, 42, 44 by performing another summation 
in the form of an 'Integral from 0 to y' along the reservoir, as follows: £ 

In figure 7, in connection with the Pressure Boundary Integrand 46, recall from figure 10 
that a 'plurality of pressure values' 30, 32, and 34 have each been calculated in ID along 
the x-axis at a 'single point in space' and at a 'single point in time' for the plurality of 
formation layers 36a, 36b, and 36c and that an 'additional plurality of pressure values' 40, 
42, 44 have each been calculated in 2D along the y-axis at a 'single point in space' and at 
a 'single point in time'. In figures 7 and 1 1, the Pressure Boundary Integrand 46 of figure 
7 will add a further additional dimension to each formation layer 36a, 36b, 36c by 
converting the 'two dimensional (2D)' calculation of pressure values 30/32/34 and 
40/42/44 for each formation layer 36a, 36b, and 36c (at the 'single point in space' and the 
'single point in time') as shown in figure 10 into a 'three dimensional (3D)' calculation of 
pressure values for each formation layer 36a, 36b, and 36c (at the 'single point in space' 
and the 'single point in time') as shown in figure 1 1 . 
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In figure 1 1, for example, the Pressure Boundary Integrand 46 of figure 7 will calculate a 
'further additional plurality of pressure values' 48, 50, 52 along the z-axis of each 
formation layer 36a, 36b, 36c, each pressure value of the 'further additional plurality of 
pressure values' 48, 50, 52 being calculated at a 'single point in space' and at a 'single 
point in time'. Therefore, in figure 1 1, the 'plurality of pressure values' 30, 32, 34 have 
each been calculated in ID in connection with the x-axis in each formation layer 36a, 
36b, 36c at a 'single point in space' and at a 'single point in time'; the 'additional 
plurality of pressure values' 40, 42, 44 have each been calculated in 2D in connection 
with the y-axis in each formation layer 36a, 36b, and 36c at a 'single point in space' and 
at a 'single point in time'; and, in addition, the 'further additional plurality of pressure 
values' 48, 50, 52 have each been calculated in 3D in connection with the z-axis in each 
formation layer 36a, 36b, and 36c at a 'single point in space' and at a 'single point in 
time'. The Pressure Boundary Integrand 46 of figure 7 actually calculates the 'further 
additional plurality of pressure values' 48, 50, 52 by performing still another summation 

in the form of an 'Integral from 0 to z' along the reservoir, as follows: £ . 

In figure 7, in connection with the Pressure Initial Integrand 54, recall from figure 1 1 that 
a 'plurality of pressure values' 30, 32, and 34 have each been calculated in ID along the 
x-axis at a 'single point in space' and at a 'single point in time' for the plurality of 
formation layers 36a, 36b, and 36c; that an 'additional plurality of pressure values' 40, 
42, 44 have each been calculated in 2D along the y-axis at a 'single point in space' and at 
a 'single point in time' for the plurality of formation layers 36a, 36b, and 36c; and that a 
'further additional plurality of pressure values' 48, 50, 52 have each been calculated in 
3D along the z-axis at a 'single point in space' and at a 'single point in time' for the 
plurality of formation layers 36a, 36b, and 36c. In figure 1 1, note that all of the pressure 
values 30 through 52 along the x, y, and z axes have been calculated 'at a single point in 
time'. The whole purpose of the 'Prediction' step 1 8 of figure 4 is to produce a 
prediction of the pressure in the reservoir (and other production data) at any point in 
space and at 'any point in time' in the reservoir. Therefore, there is a further need to 
produce a prediction of the pressure in the reservoir at 'any point in time', which includes 
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any future time. In figures 7 and 1 1, the Pressure Initial Integrand step 54 of figure 7 will 
receive the 'plurality of pressure values 30-52' for the formation layers 36a, 36b, 36c 
shown in 3D in figure 1 1 (where each pressure value exists at a 'single point in space' 
and a 'single point in time') and, responsive thereto, the Pressure Initial Integrand step 54 
of figure 7 will generate a 'time dependent plurality of pressure values' where each 
pressure value of the 'time dependent plurality of pressure values' will exist not only at a 
'single point in space' but also at 'any point in time'. Therefore, in figures 7 and 1 1, the 
'output' of the Pressure Initial Integrand step 54 will be comprised of 'one or more 
pressure values' 30, 32, 34, 40, 42, 44, 48, 50, 52 associated with 'one or more locations' 
along three dimensions in the reservoir, where: (1) 'each pressure value' of the 'one or 
more pressure values' will have a 'specific point in space' associated with the 'one or 
more locations' along three dimensions in the reservoir, and (2) 'each pressure value' of 
the 'one or more pressure values' will exist at a 'specific point in time' associated with 
the 'specific point in space' in the reservoir (that is, at any future time). The Pressure 
Initial Integrand 54 of figure 7 actually calculates the 'one or more pressure values', 
where each pressure value has a 'specific point in space' and will exist at a 'specific point 
in time', by performing still another summation in the form of an 'Integral from 0 to t' 

along the reservoir, as follows: [ . 



In figure 7, in step 56, recall the summations in the form of integrals £ £ £ Jj 

The pressures values 30-52 were calculated (in '3D' and as a function of 'time') by a 
successive integration of a Point Source Solution in the Diffusivity Equation where the 
flow rate is governed by Darcy's Law. The Gas Reservoir Evaluation and Assessment 
Tool (GREAT) of the present invention, as shown in figure 3, offers a complete workflow 
for rapid evaluation of gas reservoirs. The GREAT tool includes the Analytical Engine 20 
which further includes a newly derived solution of diffusivity equations for multiple wells, 
horizontal or vertical, in a single phase layered system under a variety of boundary 
conditions. The solution of these equations model both transient and steady state flow 
regimes and is applicable to both testing and long term performance prediction. The 
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equations applicable to laminar flow of fluids in a porous medium were the results of 
Darcy's experimental study of the flow characteristics of sand filters. This combined with 
the equation of continuity and an equation of state for slightly compressible fluid yields the 
diffusivity equation, which is the equation for pressure diffusion in a porous medium. 
Solution of the diffusivity equation under different boundary conditions forms the basis for 
prediction of bottom hole pressure response of a producing well. These analytical solutions 
are generally applicable for a single well and are used widely in the area of well testing. 
The efficiency of analytical models is generally judged by accuracy and speed. The novel 
set of solutions used in the GREAT tool of figure 3 is applicable to multiple wells, which 
can be vertical as well as horizontal. These wells can be operating as producers or injectors 
thus being of additional significance to gas well storage. The solutions have been derived 
by application of successive integral transforms. The application of these new solutions is 
characterized by stability and speed. 

Referring to figure 12, an example of the Output Record 22bl generated by the Recorder 
or Display Device 22b of figure 3 is illustrated. This Output Record 22bl provides the 
Prediction 18 illustrated in figure 4, the Output Record 22bl illustrating one or more 
pressure values (and/or other production data), each pressure value existing at any 'single 
point in space' along three dimensions (3D) in a reservoir, each pressure value existing at 
'any future point in time' associated with a 'single point in space' in the reservoir. 

A functional description of the operation of the Gas Reservoir Evaluation and Assessment 
Tool (GREAT) 22 of the present invention shown in figure 3 will be set forth in the 
following paragraphs with reference to figures 1 through 12 of the drawings. 

In figure 3, the Gas Reservoir Evaluation and Assessment software 22cl is loaded from a 
CD-Rom (or other such program storage device) into the memory or program storage 
device 22c of the Gas Reservoir Evaluation and Assessment Tool (hereinafter, the 
'GREAT tool') shown in figure 3 for storage therein. In figure 3, the Input Data 24 is 
also input to the workstation 22 of figure 3 which represents the GREAT tool of the 
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present invention. In figure 5, the Input Data 24 includes the following 'particular data' : 
well radius, skin, perforation length, fluid properties, reservoir description, and rate data. 
The aforementioned 'particular data' represents a reservoir, such as a gas reservoir, which 
includes one or more wellbores. In figure 3, the processor 22a executes the Gas 
Reservoir Evaluation and Assessment software 22c 1 stored in the memory/program 
storage device 22c and, responsive thereto, the processor 22a will generate an Output 
Record 22b 1, the Output Record 22bl being recorded or displayed on the Recorder or 
Display device 22b of figure 3. When the processor 22a executes the Gas Reservoir 
Evaluation and Assessment software 22c 1, the steps shown in figure 4 will be executed in 
sequence. In figure 4, in response to the Input Data 24, the processor 22a begins the Data 
Entry step 10 of figure 4 wherein the Input Data 24 will be received into the processor 
22a. At this point, in figure 4, the processor 22a will commence the Analytical Engine 
step 20 wherein the Analytical Engine 20 will utilize the Input Data 24 to generate a 
'computer model' 20a. In figure 4, the processor 22a will now commence the 
Verification step 14 wherein the 'computer model' 20a is verified to determine whether 
the 'computer model' will generate accurate results. In figure 6, during the Verification 
step 14, the 'computer model' 20a will undergo verification via the Test Interpretation 
step 14a and the History Matching step 14b. In both the Test Interpretation step 14a and 
the History Matching step 14b, previously known 'historical data' (having 'known 
historical results') will be introduced into the 'computer model' 20a. Responsive thereto, 
the 'computer model' 20a will generate 'results'. The 'results' generated by the 
'computer model' 20a will be compared to the 'known historical results'. If the 'results' 
approximately equal the 'known historical results', the 'computer model' 20a has 
successfully passed the verification test embodied in the Verification step 14. In figure 4, 
at this point, the processor 22a will now commence the Prediction step 18 wherein the 
future behavior of the reservoir can be predicted. During the Prediction step 18, in 
response to certain input data which may include a 'particular point in space' in the 
reservoir and a 'particular future point in time', the processor 22a will predict the pressure 
in the reservoir (as well as other production data) at 'the particular point in space' in the 
reservoir and at 'the particular future point in time'. As a result, the processor 22a can 
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predict the future behavior of the reservoir at that 'particular point in space' in the 
reservoir and at that 'particular future point in time'. The prediction that is generated by 
the processor 22a during the Prediction step 18 (i.e., the pressure and other production 
data in the reservoir at 'the particular point in space' and at 'the particular future point in 
time') is reflected on the Output Record 22bl which is recorded or displayed on the 
Recorder or Display device 22b of the GREAT tool shown in figure 3. 

Referring to figures 4 and 7 through 1 1 , recall from figure 4 that the processor 22a will 
practice the Analytical Engine step 20. During the Analytical Engine 20 step of figure 4, 
the Input Data 24 is utilized by the processor 22a to generate a 'computer model' 20a. 
The 'computer model' 20a is generated by the processor 22a during the Analytical Engine 
step 20 in the following manner. 

In figure 7, during the Analytical Engine step 20, the Input Data 24 is received for use 
during the Model Initialization step 12a. In the Model Initialization step 12a, in response 
to the Input Data 24, a pressure calculator 12a will calculate a plurality of pressure values 
within a single layer of a reservoir in one dimension (ID), such as along the x-axis, at a 
'single point in space' and at a 'single point in time'. For example, figure 8 depicts a 
reservoir wherein the plurality of pressure values 28 within a single layer of the reservoir 
has been calculated by the pressure calculator 12a in one dimension (ID), such as along 
the x-axis, at a 'single point in space' and at a 'single point in time'. In figure 7, the 
Multi Layer Loop step 26 will calculate a plurality of pressure values within multiple 
layers of a reservoir in one dimension (ID), such as along the x-axis, at a 'single point in 
space' and at a 'single point in time'. For example, figure 9 depicts a reservoir wherein 
the plurality of pressure values 30, 32, 34 within multiple layers 36a, 36b, 36c of the 
reservoir has been calculated in one dimension (ID), such as along the x-axis, at a 'single 
point in space' and at a 'single point in time'. In figure 7, the Pressure Source Integrand 
step 38 will calculate a plurality of pressure values within multiple layers of the reservoir 
in two dimensions (2D), such as along the x-axis and the y-axis, at a 'single point in 
space' and a 'single point in time'. For example, figure 10 depicts a reservoir wherein the 
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plurality of pressure values 30, 32, 34 along the x-axis of the multiple layers 36a, 36b, 
36c and the plurality of pressure values 40, 42, 44 along the y-axis of the multiple layers 
36a, 36b, 36c of the reservoir has been calculated in two dimension (2D), such as along 
the x-axis and the y-axis, at a 'single point in space' and at a 'single point in time'. In 
figure 7, the Pressure Boundary Integrand 46 will calculate a plurality of pressure values 
within multiple layers of the reservoir in three dimensions (3D), such as along the x-axis 
and the y-axis and the z-axis, at a 'single point in space' and a 'single point in time'. For 
example, figure 1 1 depicts a reservoir wherein the plurality of pressure values 30, 32, 34 
along the x-axis of the multiple layers 36a, 36b, 36c and the plurality of pressure values 
40, 42, 44 along the y-axis of the multiple layers 36a, 36b, 36c and the plurality of 
pressure values 48, 50, 52 along the z-axis of the multiple layers 36a, 36b, 36c of the 
reservoir has been calculated in three dimension (3D), such as along the x-axis and the y- 
axis and the z-axis, at a 'single point in space' and at a 'single point in time'. In figure 7, 
the Pressure Initial Integrand 54 will calculate the plurality of pressure values within 
multiple layers of the reservoir in three dimensions (3D), such as along the x-axis and the 
y-axis and the z-axis. Each of the plurality of pressure values calculated by the Pressure 
Initial Integrand 54 is associated with a 'single point in space'. However, the Pressure 
Initial Integrand 54 will also calculate the plurality of pressure values within multiple 
layers of the reservoir in three dimensions (3D) at 'any future point in time'. 

Detailed Description of the Invention 

The 'Detailed Description of the Invention' portion of this specification, set forth below, 
provides a 'detailed specification document' which discloses the construction of the 
Analytical Engine 20 of figure 4. That 'detailed specification document', which 
represents a 'compendium of analytic solutions', is set forth below, as follows. 
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Mathematical Op rations of Special Functions 

1 The diffusion mode of transference of heat, mass and pressure. 

1.1 The flux laws and the differential equation of diffusion in one dimension 

1.2 Pressure diffusion in a three dimensional porous medium 

2 Integral transforms and their inversion formulae 

2.1 Laplace Transform 

2.2 Fourier TVansforms 

2.3 Finite Fourier TVansforms. . . 

2.4 Hankel and Weber TVansforms < 

2.5 Finite Hankel TVansforms 

3 Semi-infinite medium 

3.1 The medium is bounded by the plane x = 0 and extends to oo in the direction of x positive. 
At x = 0, p(0 t t) = 0, for all t > 0 and p tends to zero as x oo. The initial pressure 
P(*,0) = 0 

3.2 The problem of 3.1, except at t = 0, p(x,0) « y>(x); <p(x) and its derivative tend to zero as 
x -» oo 

3.3 The problem of 3.1, except at t = 0, p{x t 0) = p Jt a constant for all x > 0 

3.4 The problem of 3.1, except at x = 0, p (0, t) = $ (t), an arbitrary function of time. The initial 
pressure p(x, 0) = 0 

3.5 The problem of 3.1, except at x = 0, p (0, t) = rp (*), an arbitrary function of time. The initial 
pressure p(x, 0) = (p(x)\ y>(x) and its derivative with respect to x tend to zero as x oo. 

3.6 The problem of 3.1, except at x = 0, p (0, t) = ^ ('*), an arbitrary function of time. The initial 
pressure p(x, 0) = p/ t a constant for all x > 0 

3.7 The problem of 3.1, except at x = 0, there exist a no-flow boundary, that is, ^j^2 = 0. 

-» 0 as x -+ oo. The initial pressure p(x, 0) = 0 

3.8 The problem of 3.7, except at t = 0, p(x,0) = ip{x)\ y>(x) and its derivative tend to zero as 

X 00 

3.9 The problem of 3.7, except at t = 0, p(x, 0) = p h a constant for all x > 0 

3.10 The problem of 3.7, except at x = 0, ^1 = - (£ ) ^(t), an arbitrary function of time. The 
initial pressure p(x,0) = 0 

3.11 The problem of 3.7, except at x = 0, = _ (g) ^ arbitrary f unc tion of time. The 
initial pressure p(x,0) = y?(x); <p(x) and its derivative with respect to x tend to zero as x oo. 
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3.12 The problem of 3.1, except SegH . - (g)*( t ), an arbitrary function of time. The initial 
pressure p(x, 0) = p/, a constant for all i > 0 

3.13 The problem of 3.1, except _ x p (0,t) = 0. The initial pressure p(x,0) = o! ...... 

3.14 The problem of 3.13, except at t = 0, p(x,0) = ip(x); <p(x) and its derivative tend to zero as 

X -* 00 

3.15 The problem of 3.13, except at t = 0, p(x, 0) = p/, a constant for all x > 0. ........ 

3.16 The problem of 3.13, except at x = 0, _ X p(0,t) = - (jf)^(t), an arbitrary function' 
of time. The initial pressure p(x,0) — 0 

3.17 The problem of 3.13, except at x = 0, ^gfl - Ap (0, Q = - (£) 0 (*), an arbitrary function 
of time. The initial pressure p(x, 0) = <p(x); y>(x) and its derivative with respect to x tend to 
zero as x — ► oo 

3.18 The problem of 3.1, except at x = 0, Agfl - Ap'(0, <) = - (ff)V(0/u "arbitrary function of' 
time. The initial pressure p(x, 0) = p/, a constant for all x > 0 

I Bounded medium 

4.1 The medium is bounded by the planes x = 0 and x = a. p(0,*) = p(a,*) = 0. The initial 
pressure p (x, 0) = 0 

4.2 The problem of 4.1, except at t = 0, p(x, 0) = v?(x), an arbitrary function of x, 0 < x < a. . 

4.3 The problem of 4.1, except at the extremities x = 0 and x = a, the pressures are (0 and 
ip a arbitrary functions of time. The initial pressure p(x, 0) = 0,0 < x < a 

4.4 The problem of 4.1, except at the extremities x = 0 and x = a, the pressures are ip Q (t) and 

(0, arbitrary functions of time. The initial pressure p(x,0) = v?(x), an arbitrary function 
of x, 0 < x < o 

4.5 The problem of 4.1, except at x = 0, p(0, t) = 0 and at x = a, = 0 . The initial pressure 
p(x, 0) = 0, 0 < x < a. 

4.6 The problem of 4.5, except at x = 0, p(0, 0 = 0 and at x = a, = 0 . The initial pressure 
p(x,0) = v?(x), an arbitrary function of x,0 < x < a 

4.7 The problem of 4.5, except at x = 0, p(0, t) = tfo(0 and at x = o, = - (f ) 0 a (t). (0 
and ^ a (£) are arbitrary functions of time. The initial pressure p(x,0) = 0, 0 < x < a 

4.8 The problem of 4.5, except at x = 0, p(0, t) = ^ 0 {t) and at x = o, = - (£) ^o(0- (0 
and Vo(0 are arbitrary functions of time. The initial pressure p(x, X 0) = y?(x), an arbitrary 
function of x,0 < x < a 

4.9 The problem of 4.1, except at x = 0, p(0, t) = 0 and at x = a, + A a p (a, 0 = 0. The 
initial pressure p(x, 0) = 0, 0 < x < a 

4.10 The problem of 4.9, except at x = 0, p(0,t) = 0 and at x = a, -f A a p(a,*) = 0. The 
initial pressure p(x, 0) = y?(x), an arbitrary function of x, 0 < x < a? 

4.11 The problem of 4.9, except at x = 0, p(0, t) = ^ 0 (0 and at x = o, 

^^H* + A a p (a, t) = - (£ ) 0 a (t). (0 and (i) are arbitrary functions of time. The initial 
pressure p (x, 0) = 0, 0 < x < a 

4.12 The problem of 4.9, except at x = 0, p(0,*) = tp Q (t) and at x = a, 

+ A a p (a, t) = - (£) tf a (*)• (0 and r/> a (t) are arbitrary functions of time. The initial 
pressure p (x, 0) = tp(x)j an arbitrary function of x, 0 < x < a 

4.13 The problem of 4.1, except at the extremities we have no flow boundary conditions, that is 
at x = 0, -3&a = 0 and at x = a, = 0. The initial pressure p(x,0) =0, 0 < x < a. /. 

4.14 The problem of 4.13, except at x = 0, ^ = 0 and at x = a, = 0. The initial 
pressure p(x, 0) = y?(x), an arbitrary function of x, 0 < x < o 

4.15 The problem of 4.13, except at x = 0, ^1 = - (g) ^ Q and at x = a, 

^f 1 ^ = ~ (* ) tfa (0- W and (0 are arbitrary functions of time. The initial pressure 
P(*,0) = 0, 0<x<o 
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416 S:, t r blem ° f 413 ' mS * at X = °' = - (f ) * W and at x = a, 

X ~ ^ a W\ ^° ( *> 40(1 W *"> arbitrary functions of time. The initial pressure 
p(x, 0) = <p(x), an arbitrary function of x,0 < x < a T h 

417 J(t!o p ) ro = bl o em o< i < T pt °' * = 0 and at 1 = °- p {a ' 0 - a The initial ' 

4.18 The problem of 4.17 except at t = 0, p(x,0) = ^(x)/an arbitrary* function of x,0 < x < a' 

4.19 The problem of 4.17, except at x = 0, *g& = _ (f j ^ (t) ^ at x ( t) = ( " " 
1>o (*) and xl> a (t) are arbitrary functions of time. The initial pressure p(x,0) = 0 0 < x < o 

4.20 The problem of 4.17, except at x = 0, ^ = - (f ) ^ (,) and at x = a, p (a', t) = tf a (t) ' 
Vo i t) and V« (t) are arbitrary functions of time. The initial pressure p(x,0) = v(x) an 
arbitrary function of x, 0 < x < a ^ ; ' 

4.21 The problem of 4.1, except at x = 0, = 0 and at x = a Sei&A+\' n(a t \l n " T v " 
initial pressure p(x, 0) = 0, 0 < x < «.7T ... * ( ' " 

4.22 The problem of 4.21, except at x = 0, flgfl = 0 and at x'= a, + A 0 p(o,t) = 0 The' 
initial pressure p(x, 0) = <p{x), an arbitra? function of x, 0 < x < o 

4-23 SIS' 0 " 6111 ° f 4 ' 21 ' eXCCPt at 1 = °> ^ = " (f ) (0 and at x = a, 

+ X f fe' ^ = " (P ^ a W> ^° (<) and ^ a W are arbitrary functions of time. The initial 
pressure p(x, 0) = 0, 0 < x < o 

4 24 SS.0 rOblem ° f 4 ' 21 ' CXCePt at x = °> ^ = - (f ) ^ (t) and at x = o, 

+ A /° P fc' 4) = , ^ a ^° W and ^° W we arbitrary functions of time. The initial 
pressure p(x, 0) = ^(x), an arbitrary function of x, 0 < x < a 

4.25 The Problem of 4.1, except at x = 0, SE^Q. - A o p(0, *) = 0 and at x = a, 

+ A a p (a, = 0. The initial pressure p(x, 0) = 0, 0 < x < a 

4 26 Slo r ° blem ° f 4 ' 25 ' at 1 = °- ^ " A 0P (0, «) = 0 and at x = a, 

+ A a p (a, t) = 0. The initial pressure p(x, 0) = ip(x), an arbitrary function of x 
0 < x < a ' 

4 27 2ST 0blem ° f 4 ' 25 ' 8X0ept at 1 = °' ^ - A 0P(0, t) = - (jj) (f j and at x = 'a, 

+ Vfc f) = " (P ^ a ( ^ ^° (<) and ^ arbitrary functions of time. The initial 
pressure p(x, 0) = 0, an arbitrary function of x, 0 < x < a 

4 ' 28 SlS r ° blem ° f 4-25 ' eXCCPt &t X = °' ^ " AoP (°- *> = - («) *> (0 and at x = a, 

+ V ^ = ,~J^ ^° ^° W and Va(<) are arbitrary functions of time. The initial 
pressure p(x, 0) = y>(x), an arbitrary function of x, 0 < x < a 

4.29 The problem of 4.1 except at x = 0, _ Ao p( 0) t) = 0 and at x = a ',p(a,t) = 0. The' 
initial pressure p(x,0) = 0, 0 < x < a 

4.30 The problem of 4.29 except at x = 0, ^ _ Aop((U ) . 0 a nd at x = a, p(a, i) - 0. The' 
initial pressure p(x, 0) = v (x), an arbitrary function of x, 0 < x < a 

4.31 The problem of 4.2S », except at x = 0, ^ - A 0 p (0, <) = - (f) (t) and at x ='a, 
p o, tj = rp a {t). i> 0 (t) and Va(r) are arbitrary functions of time. The initial pressure 
p(x,0) = 0, 0<x<4 v 

4.32 T fa e Problem of 4.29 except^ 

/«%*■/ ?■ ^° *«W are arbitrary functions of time. The initial pressure 
p(x,0) = <p(x), an arbitrary function of x,0 < x < o 

4.33 The problem of 4.1, except at x = 0, - A 0 p(0,t) = 0 and at x = a, ^ = 0 The 
initial pressure p(x,0) = 0, 0 < x < a °* 

4.34 The problem of 4.33, except at x = 0, ^ - A o p(0,*) = 0 and at x = o, 2« = 0 The 
initial pressure p(x, 0) = ^(x), an arbitrary function of x, 0 < x < a .... 
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4.35 ^Problem of 4.33, except at * « 0, ^ - A oP (0, — (0 and at * =* a, 

*5> : 0 : 2&*iS2£^l,r.?^ ta-taB of time - The ^ ^ 

4.36 ^Problem of 4.33 except at x = 0, _ Aop (0 , t) = _ ^ ^ (t) ^ x ^ 

irfwrt ~ ( t W * aad ^ K^tty Unctions of time. The initial pressure 
p(x, 0) = v>(x), an arbitrary function of x, 0 < x < a pressure 

5 Cascaded medium 

5.1 N connected mediums a, < x < a j+u V j = 0, 1 , K — 1. At x = ao 

= -(«*<*> and at x = a*, ^ = _(jf)* (t) . At the interface x = 0j) 
= 1,2 * " l > * W = " ft), = " (^#^) and 

c o 2^ W J7 ^ Pi " 1 ' ^ " p > (o >'' *)},*> 0. The initial pressure p, (x, 0) = 0 

5.2 The problem of 5.1 except, the initial pressure Pj (x, 0) = ^(x) 

5.3 The problem of 5.1 except, at x = a«, p (a„, t) = ^ (t). The initial pressure p/(x, 0) - ^(x)' 

54 Sls^^t^ 1 : r.r - ***-^ ( -"^ - 

5 5 SLSm,, 0) t 1 ^. "T^. that ^ . thC lMt mCdiUm W Semi '- infinite - ' T «« initial • 

5.6 The problem of 5.1 exrent at t _ „ i\ /. /a' . ' 

and at x = an, 




5. i--— «• ~«fv, ai x = ao, p lfl0 , tj = Vo (t) and at x = on, p (o K , t) = tfo ft) The 
initial pressure pj- (x,0) =vjj(x) 1 w 

5-8 SlS 0blem ° f 51 eXCept> at 1 = °o. P («o. *) = W and at x = a«, 

c « t X " P ( °* ' l) = " The initial Pressure p.- (x, 0) = u>Ax). . 

ISfn at 1 = 001 P / ao ' 0 = ^ W and °« ^ °°- tha < ^, the last medium ' 

is semi-infinite. The initial pressure pj (x,0) = <pj(x) 

5.10 Theproblem of 5.1 except, at x = a 0 , _ XqJ> {0>t) m 'J (f j ^ " ^ ^ , = ^ ' 

e* = ~ ( %) V>o„ (<)• The initial pressure Pj (x, 0) = ^(x) 

5.11 The Problem* 51 except at x = a 0 , ^ -A 0 p(0,<) = - (f)^( t )and at x'= o N ' 
P(a H ,t) = tl^ (t). The initial pressure Pj (x, 0) = ^ (x) 

512 SLStr ^ 6XC t\ at X = a °' ^ - A ° P(M) = -(V)^W and at . - a N ,' 

ei + A N p (o n , t) = - (£) ^, OK (<). The initial pressure Pj (x, 0) = <pj(x) 

513 ^I^H- ° f 51 eXCePt, / t X = ?' ~ AoP (0 ' *) — (!)* W - oo," that is, ' 
the last medium is semi-infinite. The initial pressure Pj (x,0) = ^(x). .......... .7 

Semi-infinite quadrant 

61 Jf^ n n l iUm " J° Unded , by th A e P lanes « - 0 and » - 0j « and y extend to oo in the directions 
of x pos,t,ve and y pos.t.ve. At x = 0, p (0, y, t) = 0, y > 0, t > 0 and at „ = 0, p (x, 0, t) = 0 

fi 2 n 'i ^ t °" , °" , " 00Md »- >fla The initial pressure p(x,yb) = 0 . 

JnH, P ; ,T CXCePt ' P(l,y,0) = and its derivativ es with respect to x 

and y tend to zero as x -* oo and y -» oo respectively 

6.3 The problem of 6.1 except, p(x, y, 0) = P/ a constant for all x > 0 and y > 0 

6.4 The problem of 6.1, except for all t > 0, p (0, y, t) = ^ y (y, t ), y > 0 and p (x, 0, 1) = ^ (x,' t), ' 

and ' T PrCSSUre P(X ' I' 0) = V) ' y) " d itS derivatives respect to . 
and y tend to zero as x ► oo and y — > oo respectively. . 

6.5 The problem of 6.4, except, p(x, y, 0) = P/ a constant for all x > 0 'and y > 0 

6.6 The problem of 6.1, except, = 0 and p(0,y,0 = 0. The initial pressure is zero. .' .' .' .' 
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6.7 The problem of 6.6, except, p(x, y, 0) = <p(x, y); y>(x, y) and its derivatives with respect to x 
and y tend to zero as x oo and y — ► oo respectively. 

6.8 The problem of 6.6, except, p(x, y, 0) = p, a constant for all x > 0 and y > 0. . ! 

6.9 The problem of 6.6, except for all t > 0, = _ ^ j ^ x ( x , t ), x > 0 and p (0, y, t) = 

V'vtv.O. V > 0. The initial pressure p(x,y,0) = <p(x,y); y>(^, y) and its derivatives with 
respect to x and y tend to zero as x oo and y — ► oo respectively. 

6.10 The problem of 6.9, except, p(x, y, 0) = p/ a constant for all x > 0 and y > 0 

6.11 The problem of 6.1, except, ^i^i _ Ap(x,0,*) = 0 and p(0,y,i) = 0. The initial pressure 
is zero 

6.12 The problem of 6.11, except, p(x,y,0) = <p(x,y); v?(x,y) and its derivatives with respect to x 
and y tend to zero as x — ► oo and y — ► oo respectively. . . . 

6.13 The problem of 6.11, except, p(x, y, 0) = p/ a constant for all x > 0 and y > 0. 

6.14 The problem of 6.11, except for all t > 0, _ Xp ^ 0| t) - - (^) ^ x (x, t) f x > 0 and 
P (0, y, t) = ^ v (y, t), y > 0. The initial pressure p(x, y, 0) = y?(x, y); <p(x, y) and its derivatives 
with respect to x and y tend to zero as x oo and y — > oo respectively. 

6.15 The problem of 6.11, except at t = 0, p(x,y,0) = p/ a constant for all x > 0 and y > 0 

6.16 The problem of 6.1, except, at x = 0, 3^*1 = 0 and at y = 0, = Q The MfM 
pressure is zero. 

6.17 The problem of 6.16, except, p(x,y,0) = <?(x,y); v?(x,y) and its derivatives with respect to x 
and y tend to zero as x — ♦ oo and y — > oo respectively. 

6.18 The problem of 6.16, except, p(x, y, 0) = p 7 a constant for all x > 0 and y > 0 

6.19 The problem of 6.16, except for all t > 0, = _ ^ ^ ( y>t ) j y > 0 and ^E^M = 

~ (£) iM*.*). * > 0. The initial pressure p(x,y,6) = y>(x,y); y?(x,y) and its derivatives 
with respect to x and y tend to zero as x — ► oo and y -* oo respectively. 

6.20 The problem of 6.16, except, p(x, y, 0) = p/ a constant for all x > 0 and y > 0 

6.21 The problem of 6.1, except, 2e^1 = 0 and at y = 0, - Ap (x, 0, t) = 0. The initial 
pressure is zero 

6.22 The problem of 6.23, except, p(x, y, 0) = <p(x, y); y?(x, y) and its derivatives with respect to x 
and y tend to zero as x — ► oo and y — ► oo respectively. 

6.23 The problem of 6.23, except, p(x, y, 0) = p 7 a constant for all x > 0 and y > 0 

6.24 The problem of 6.21, except for all t > 0, = - (^) ifi y (y,*), y > 0 and _ 

Ap(x,0,t) = -(^-) ipx(x>t), x > 0. The initial pressure p(x,y,0) = y?(x,y); <p(x,y) and its 
derivatives with respect to x and y tend to zero as x -> oo and y oo respectively. 

6.25 The problem of 6.24, except at t = 0, p(x, y, 0) = p/ a constant for all x > 0 and y > 0 

6.26 The problem of 6.1, except, at x = 0, - A x p(0,y,*) = 0 and at y = 0, _ 
A y p(x,0,£) = 0. The initial pressure is zero 

6.27 The problem of 6.26, except, p(x,y,0) = <p{x y y)\ v?(x,y) and its derivatives with respect to x 
and y tend to zero as x — ► oo and y — ► oo respectively. 

6.28 The problem of 6.26, except, p(x, y, 0) = p/ a constant for all x > 0 and y > 0 

6.29 The problem of 6.26, except for all t > 0, - A x p(0,y,*) = - tf»(y,t), V > 0 
and gE^M _ A y p(x,0,*) = - tf> x (x,<), x > 0. The initial pressure p(x,y,0) = <p(x,y); 
<p(x, y) and its derivatives with respect to x and y tend to zero as x oo and y oo 
respectively. 

6.30 The problem of 6.26, except, p(x, y, 0) = p/ a constant for all x > 0 and y > 0 
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7 Semi-infinite lamella , 

7.1 The medium is bounded by the planes x = 0, y = 0 and y = 6; x — oo in the direction of x 
positive. p(0,y,t) = 0, for 0 < y < 6, t > 0 and p(x,0,t) = p(x,M) = 0, for all x > 0, t > 0. 
p tend to zero as x — ► oo. The initial pressure is zero 

7.2 The problem of 7.1, except at t = 0, p(x, y, 0) = <p(x t y); <p(x, y) and its derivative with respect 
to x tend to zero as x — ► oo 

7.3 The problem of 7.2, except, p(x, y, 0) = p 7 ; p 7 , a constant for all x > 0, 0 < y < 6 

7.4 The problem of 7.1, except, t > 0, 0 < y < 6, p(0,y,*) = V y (y,*) and x > 0, p(x,0,0 = 
tfeoCM). P(x,M) = The initial pressure p(x,y,0) = <p{x,y)\ y?(x,y) and its 
derivative with respect to x tend to zero as x — > oo 

7.5 The problem of 7.4, except at t = 0, p(x, y, 0) = p /; p 7 , a constant for all x > 0, 0 < y < 6. . . 

7.6 The problem of 7.1, except, 2el*M = _ ^ ^ (Xj t)| p (Xj 6> t) = ^ (x> ^ ^ 

p(0,y,*) = tl> y {y t t). The initial pressure p(x,y,0) » y>(x,y); <p(x,y) and its derivative with 
respect to x tend to zero as x — + oo 

7.7 The problem of 7.6, except, the initial pressure p (x, y, 0) = p/, a constant 

7.8 The problem of 7.1, except, _ Ax0p (a5| 0 , i) = - (^) V*o (x, *), P (x, 6, 0 = 4.6 (x, 0. 

and p(0,y,*) = V>y (y,0- The initial pressure p (x, y, 0) = <?(x,y); <p(x,y) and its derivative 
with respect to x tend to zero as x — > oo 

7.9 The problem of 7.8, except, p(x, y, 0) = p 7 , a constant 

7.10 The problem of 7.1, except, p (x, 0, t) = 0 xO (x, 0, = - (#) 4*6 (x, *), and p (0, y, 0 = 

1> y (y,i). The initial pressure p(x,y,0) = y?(x,y); y?(x,y) and its derivative with respect to x 
tend to zero as x — > oo 

7.11 The problem of 7.10, except, p(x, y, 0) = p /f a constant 

7.12 The problem of 7.1, except, = _ (^)^ x0 (x,0, = - (£)lMM) and 

P(0,y,t) = V v (y,0- The initial pressure p(x,y,0) = y?(x,y); y?(x,y) and its' derivative with 
respect to x tend to zero as x — ♦ oo 

7.13 The problem of 7.12, except p(x,y,0) = p 7 , a constant 

7.14 The problem of 7.1, except, *^- AxO p(x,0, <)sB5 _ ^ (Mf SegUU. 

and p (0, y, t) =j> v (y, t). The initial pressurep(x, y, 0) = tp(x t y); y>(x, y) and its derivative with 
respect to x tend to zero as x — ► oo 

7.15 The problem of 7.14, except p(x, y, 0) = p 7 , a constant 

7.16 The problem of 7.1, except, p (x, 0, t) = V*o (x, 0, 5e %^ + A* 6 p (x, 6, <) = - ^ x6 (x, £), 

and p(0,y,*) = il> y (y t t). The initial pressure p(x, y, 0) = y?(*.y); v>(x,y) and its derivative 
with respect to x tend to zero as x — ► oo 

7.17 The problem of 7.16, except p(x,y,0) = p/, a constant 

7.18 The problem of 7.1, except, = _ ^ ^ ^ 

9j> % b,t} + <W(x,M) = - (£)^ x &(x,0 and p(0,y,i) = 4 v (y,*)- The initial pressure 
p(x, y, 0) = y?(x, y); v?(x, y) and its derivative with respect to x tend to zero as x oo 

7.19 The problem of 7.18, except p(x, y, 0) = p/, a constant 

7.20 The problem of 7.1, except, _ A xO p(x,0,<) = - (^) i/> x0 (x,t), 

+ A l6 p(x,M) = - (£)^x6(x,£) and p(0,y,t) = Vyfo,*)- The initial pressure 
P(x, y, 0) = v?(x, y); y?(x, y) and its derivative with respect to x tend to zero as x oo 

7.21 The problem of 7.20, except p(x, y, 0) = p/, a constant 



30 



7.22 The problem of 7. 1, except, = - ($) A (y, t), p (x, 0, t) = (x, t) and p (x, 6, «) = 

drt (x, «)• The initial pressure p(x, y, 0) = <p(x, y); <p(x, y) and its derivative with respect to x 
tend to zero as x — ♦ oo 

7.23 The problem of 7.22, except p(x,y,0)=p/, a constant 

7.24 The problem of 7. 1, except, = - (%) VxO (x, t), p (x, 6, t) = * xb (x, t), and ^1 = 

— (j^") (V. *)• The initial pressure p(x,y,0) = <p{x,y); <p(x,y) and its derivative with 
respect to x tend to zero as x — ► oo 

7.25 The problem of 7.24, except p(i,y,0)= P/ , a constant 

7.26 The problem of 7. 1, except, Se^*1 _ XxQp (x> 0< t) = _ ^ (>< f)> 

p (x, 6, t) = Vxi (x, t), and 5e^*1 = _ (jl) ^ fa t y The initial pressure p(l) ^ Q) = ^ y) . 
^(x,y) and.its derivative with respect to x tend to zero as x -» oo 

7.27 The problem of 7.26, except p(x, y, 0) = p/, a constant 

7.28 The problem of 7.1, except, p(x,0,t) = (x,t), = - (^) rj, xb ( X)t ), and ^E^i = 

~(£)lMv.O- The initial pressure p(x,y,0) = y>(x,y); ^(x.y) and its derivative with 
respect to x tend to zero as x — » oo 

7.29 The problem of 7.28, except p(x, y, 0) = p/, a constant 

7.30 The problem of 7.1, except, ftfgM = - (^)^ M^l ^ - (^j ^ (x,0, and ' 

aP ai V, ° = " (V» *)• The initial pressure p(x, y, 0) = ^(x, y); y?(x, y) and its derivative 
with respect to x tend to zero as x — > oo 

7.31 The problem of 7.30, except p(x,y, 0) = p 7> a constant 

7.32 The problem of 7.1, except, *^ - A x0 p (*, 0 , t) = - (^) ^ (*, 0, 

= - (#) 1^ (x, 0, and = _ ^ y (f/| ^ The initial presgure p(Xf y> Q) = 
^(^y); ¥>(x>y) and its derivative with respect to x tend to zero as x -* oo 

7.33 The problem of 7.32, except p(x,y,0) = p 7 , a constant 

7.34 The problem of 7.1, except, p (x, 0, t) = V*o *). 

H- A x6 p (x, 6, 0 = - [%) ^ x6 (x, t), and = _ % (|f| t y The initial pres . 

sure p(x, y, 0) = tp(x t y); y?(x, y) and its derivative with respect to x tend to zero as x oo. 

7.35 The problem of 7.34, except p(x, y, 0) = p/, a constant 

7.36 The problem of 7.1, except, = - 0 lO 

+ A x6 p (x, M) = - (£) 0 x6 ( X , t), and ft^l = _ ( jl) % (y> t)> The initial pres . 
sure p(x, y,0) = <p(x,y); y?(x,y) and its derivative with respect to x tend to zero as x -> oo. 

7.37 The problem of 7.36, except p(x, y, 0) = p/, a constant 

7.38 The problem of 7.1, except, ^^--WfeO.*)--^)^^.*), 

^ + A x6 p (x, 6, t) = - ($) ^ (x, 0, and Se^aA = (^) ^ y t) . T he initial pres- 
sure p(x, y, 0) = ip{x % y); y?(x, y) and its derivative with respect to x tend to zero as x oo. 

7.39 The problem of 7.38, except p(x,y,0) =p/, a constant 

7.40 The problem of 7.1, except, p (x, 0, t) = 0 xO *), p (x, 6, t) = V*& (*, «), and 

dp{ &z ll) - V (0» V>t) = " (^) (y, 0- Th e initial pressure p(x, y, 0) = v?(x, y); p(x, y) and 
its derivative with respect to x tend to zero as x -> oo 

7.41 The problem of 7.40, except p(x,y,0) =p 7 , a constant 
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7.42 The problem of 7. 1, except, ^E^M = _ ^ ^ (l> t)> p {x> 6 _ t) = ^ (l> ^ ^ >fr ^, t) _ 

V (0. y. 0 = " (£) 0» (y, <)• The initial pressure p(x, y, 0) = p(x, y); y) and its deriva- 
tive with respect to x tend to zero as x oo 

7.43 The problem of 7.42, except p(x,y,0) =p/, a constant 

7.44 The problem of 7.1, except, ^^i- W(*,0,t) = -(^).^«o (*,«), 

p(x,6,t) = iMM). and _ A v p(0,y,r) = -(^)^(y.t). The initial pressure 
P(x, y, 0) = y?(x, y); y?(x, y) and its derivative with respect to x tend to zero as x -» oo 

7.45 The problem of 7.44, except p(x,y,0)=p/, a constant ' ' ' 

7.46 The problem of 7.1, except, p (x, 0, t) = V«o (x, t), Se^M = _ ^ ^ (li t)> ^ ggy _ 

V (0, y, 0 = - (£) ^ y (y, t). The initial pressure p(x, y, 0) = <p(x, y); ^(x, y) and its deriva- 
tive with respect to x tend to zero as x — ♦ oo 

7.47 The problem of 7.46, except p(x, y, 0) = p/, a constant 

7.48 The problem of 7.1, except, = _ (^) ^ (a , t)| fl^a = _ ^ j ^ ^ ^ ^ " 

- A„p (0, y, t) = - ) V v (y, t). The initial pressure p(x, y, 0) = v?(x, y); y>(x, y) and 
its derivative with respect to x tend to zero as x — ► oo 

7.49 The problem of 7.48, except p(x,y,0) = p 7 , a constant 

7.50 The problem of 7.1, except, _ XxQ p (x, 0, t) = - (^) (x, *), 

= " (#) ( x > '). and ^5^1 _ Ayp (0> y> t) « _ ^ ^ (y> 0 The . n . tial pressure 
p(x, y, 0) = y>(x, y); y?(x, y) and its derivative with respect to x tend to zero as x -> oo 

7.51 The problem of 7.50, except p(x, y, 0) = p/, a constant 

7.52 The problem of 7.1, except, p(x, 0,t) = tf*o (x,*). 

^ + A x6 p (x, M) = - (£) tf» (x, t), and ^ _ AyP (0 , y% t) = _ ^ y (yj t) . The 

initial pressure p(x, y, 0) = <p(x t y); y?(x, y) and its derivative with respect to x tend to zero as 
x -* oo 

7.53 The problem of 7.52, except p(x,y,0) =p/, a constant 

7.54 The problem of 7.1, except, = - (^) ^ x0 (*, *), 

^ + A l6 p (x, M) - - ( £ ) (x, 0, and M°**2 _ XyJ} «,, y> *) = _ ( ft ) ^ {y> t) . The 

initial pressure p(x, y, 0) = y?(x, y) ; ^(x, y) and its derivative with respect to x tend to zero as 
x — ♦ oo 

7.55 The problem of 7.54, except p(x, y, 0) = p/, a constant 

7.56 The problem of 7.1, except, _ XxQ p (x, 0, t) = - (^) («, 0. 

*S*a + Ax6P (x, 6, t) = - (^) 4> xb (x, «), and ^ - x y p (0, (^) ^ (y, t). The 

initial pressure p(x, y, 0) = <p(x, y); (^(x, y) and its derivative with respect to x tend to zero as 
x — > oo 

7.57 The problem of 7.56, except p(x, y, 0) = p/, a constant 

Rectangle 

8.1 The medium is bounded by the planes x = 0, x = a, y = 0 and y = 6. p (0, y, t) = 0 and 
P(a,y,t) = 0 for 0 < y < 6. p(x,0,*) = 0 and p(x t b t t) = 0 for 0 < x < a, i > 0.' The initial 
pressure is zero 

8.2 The problem of 8.1, except at t = 0, p(x, y, 0) = <p(x, y); <p(x, y) an arbitary function of x and 
y, 0 < x < a and 0 < y < 6 
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8.3 The medium is bounded by the planes x = 0, x = o, y = 0 and y = b. p (0, y, t) •= ^ (y t) 
and p(a,y,t) = tf v „(y,t) for 0 < y < 6. p(x,0,t) = V*o(x,t) and p(x,M) = V«6(x,t) for 
0 < x < a, * > 0. The initial pressure p(x, y, 0) = <p(x, y) 1 

8.4 The problem of 8.3, except, p(0,y,*) = V v o(y,0, = - (^)i> ya (y,t), p(x,0,t) = 
1>x0 (x, t), and p (x, 6, t) = V l6 (x, t). The initial pressure p(x, y, 0) = y>(x, y) 

8.5 The problem of 8.3, except, p (0, y, t) = ^ (y, t), 2^ + A ya p (a, y, t) = - (£-) V> ya (y, t), 

p (x, 0, i) = tf> l0 (x, t) and p (x, 6, t) = (x, t). The initial pressure p(x, y , 0) = yj*(x. y) ' 

8.6 The problem of 8.3, except, p(0,y,i) = ^ 0 (y,t). - - (%) (y,<). p(x,0,*) = 
Vxo (x, t), and ^2^2 = _ ( jl) ^ ( X| t) . The initial pressure ^ „ ( 0) = ^ y) 

8.7 The problem of 8.3, except, p (0, y, = ^ y0 (y, i), = - (^) V» vo (y, t), P (x, 0, t) = 
V*o (x, *) and +A l6 p (x, *) = - ( x , i). The initial pressure p(x, y, 0) = y>(x, y). 

8.8 The problem of 8.3, except, p (0, y, f) = V y0 (y, 0, + A yo p (a, y, <) = - (^) V» yo (y, t), 

P (x, 0, t) = Vxo (x, 0 and + Art p (x, M) = - (#) ^ x6 (x, *)• The initial pressure 
P(x,y,0) = <p(x,y) " 

8.9 The problem of 8.3, except, p (0, y, t) = (y, t), p(o, y, t) = tf yo (y, f), = _ U_\ ^ (S| f)> 
and p (x, 6, i) = Vx& (x, t). The initial pressure p(x, y, 0) = yj(x, y) 

8.10 The problem of 8.3, except, p (0, y, t) = Wo (v, *). 22 ^ £t£i = - (%) 1> ya (y, «). 

ap( ly 0,t) = ~ (*;) (x, t). and p (x, b, t) = tf> xb (x, t). The initial pressure p(x, y, 0) = v?(x, y). 

8.11 The problem of 8.3, except, p (0, y, t) = ^ y0 (y, «),' + A ya p (a, y, t) = - (£-) V» yo (y, t), 
ap ^ 0,t) = ~ (*;) Vfeo (x, 0. and p (x, 6, t) = rp xb (x, i). The initial pressure p(x, y, 0) = y>(x, y). 

8.12 The problem of 8.3, except, p (0, y, *) = V„o (y, 0. P(a, y, «) = tf yo (y, <), = - ^ x0 (x, t), 
and ep( ^ 6,t) = - (£j %l> xb (x, t). The initial pressure p(x, y, 0) = <p(x, y) 

8.13 The problem of 8.3, except, p(0,y,f) = V y o(y,t). = - (fc) ^ ya (y,t), = 
" (*;) V-xO (x.t), and = - (^) Vxfc(x,i). The initial pressure p(x,y,0) = <p(x,y). . . 

8.14 The problem of 8.3, except, p (0, y, <) = V y0 (y, t), + A ya p (a, y, t) = - (^) ^ yo (y, *), 

= " (#) (x, t), and = - (^) ^ x6 («, <). The initial pressure" p(x, y, 0) = 
f(x,y) 

8.15 The problem of 8.3, except, p(0,y,t) = Mv,t),p(a,y,t) = V„„ (y,t), = _ ^(x,t), 
and + A l6 p (x, b,t) = - ) tf>*6 (x, <)• The initial pressure p(x, y, 0) = <p(x, y) 

8.16 The problem of 8.3, except, p(0,y,f) = 0 yO (y,<), = - (^)^«(y,t), gE ^ Ji = 
- (i^) V-xo (x, t), and «E^M +Aitp ( x> b> q = _ ^ ^ (i> f j The mitia , pressure p(;E> y> Q) = 
V(x,y) " 

8.17 The problem of 8.3, except, p (0, y, i) = ^ y0 (y, <). + A yo p (o, y,t) = - (^) ,p ya (y, t), 

= " (£) V-xO (x, i), and SE^i + Aitp (a:> bit) = _ ^ (I> t) The ' initial pre s. 
surep(x,y,0) = V >(x,y). . * 
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818 £S blem of 8 ' 3, except ' p (0 ' y ' 0 = {y< tl p (a> y ' l) = ( f> *>• 

- AxoP (x, 0. 0 = - (^) tfxo (x, t). and p (x, 6, <) = rp xb (y, *). The initial pressure 
P(x,y,0) = v(x,y) 

8.19 The problem of 8.3, except, p (0, y, t) = ^ (y, <), . - ^ ya ( y , t ), 

~ A *0P (*. 0, t) - - (^) Vxo (x, 0, and p (x, 6, *) = 0 l6 ( y , *). The initial pressure 
P(x,y,0) = <p(x,y) 

8.20 Theproblemof8.3,except,p(0,y,0 = Vyo(y,0, +\aP(a,y,t) = - (fij^fat), 

5E ^T ii " -W (x, 0, t) = - (£) tfxo (x, t), and p (x, 6, t) = ^6 (y, 0- The initial pressure 
p(x,y,0) = v(x,y) 

8.21 The problem of 8.3, except, p(0,y,0 = ^o(y,t).P (a, y,i) = ^o(y, t)', 

p"S" A, ° P(l,M) = "(^)^ o(x,t)l and ^ - -W *•(«.«). The initial 
P(x,y,0) = <p{x,y) 

8.22 The problem of 8.3, except, p (0, y, t) = ^ (y, «), = - ^ („, «), 

--W(x,0,t) = -(#)iM».*). and .fla^!l = -(VWx,«). The initial 
pressure ' 

P(x,y,0) =y>(x,y). . . 

8.23 The problem of 8.3, except, p (0, y, t) = rp y0 (y,t), 2e^2 + A ya p (a, y, t) = - (f)^ {y , t ), 

pS" Al ° P(X,0,<) = " &nd ^ - -(^)^6(x,0/ T he initial 
P(x,y,0) = v?(x,y) 

8.24 The problem of 8.3, except, p(0,y,t) = ^o(y,<).P (a, y,r) = ^ v « (V. 0. 

^ - W(*,0,«) - -(fc)iM*.0, and A$M + W(*,ft,t) = -f^Wx,*). 
The initial pressure p(x, y, 0) = <p(x, y) 

8.25 The problem of 8.3, except, p (0, y, t) = ^ 0 (y, t), = - (#) V-„a (y, t), 

^ - W(*,0,t) = -($)*o(»,l), and °*§M + X xb p(x,b,t) = - (*.«)■ 
The initial pressure p(x, y, 0) = y>(x, y) 

8.26 The problem of 8.3, except, p (0, y, t) = ^ (y, i), + X ya p (a, y,t) = - (g) \p ya ( y ,' t j , ' 

- AxoP(x,0,t) = - (£) and ^E^Mi + X xb p(x,b,t) - - (f) 

The initial pressure p(x, y, 0) = y>(x, y) 

8-27 The problem of 8.3, except,^ ......... 

^S 0 ' 0 = " fe) (*•'). andp(x,6,t) = i> xb (y,t). The initial pressure p(x,y,0) = y>(x,y). 

8.28 The problem of 8.3, except, = - ^ (y , 2^ — (f)*. (y, g, 

aP( ^°' <) = " (#) ^o-( x - *). 411(1 P <) = V*(> (y, <)• The initial pressure p(x, y, 0) = ^(x, y). 

8.29 Theproblemof8.3, except, = _ (^)^( y , 0i ftfcitf = _ (fc)^,,,, 
8P( ^ 0,t) = " (*;) ^xo (x, t), and p (x, 6, t) = rl> xb (y, t). The initial pressure p(x, y, 0) = <p(x, y). 

8.30 The problem of 8.3, except, = - (^) ^ ( y , ^ = _ ^ (y> t)> 

= " (#) (x, t), and £1^2 = _ fg\ ^ (li () , The initia] pressure y> Q) = 
<p(x,y) 
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8.31 The problem of 8.3, except, = - (f.) ^ ( y , t), &&*+\ yo p (o, ,,«)--. (fc) y, yo (y, t), 

= ~ (t ) <*• «). ^ = - (£) V»x6 (x, t)- The initial pressure p(x, y, 0) = 
<p(x.v) , 

8.32 The problem of 8.3, except, = - ($) ^ (y, t), &§*£+X ya p (a, y, t) = - (fc) V- vo (y, t), 

2 ^ = -(^)^xo(*,0.and52^+A lt p(« I 6, t ) = -(^)^ lt (x > 0. The initial pres- 
surep(x,y,0) =y?(x,y) - 

8.33 The problem of 8.3, except, = _ (g) ^ (y> t)i p (a _ y> t) . ^ (y> f)> 

5E ^r j ~ X *°P ( x - °« ') ™ - (#) V»iO (x, t), and p (x, 6, t) = V» l6 (y, r). The initial pressure 
p(x,y,0) = v>(x,y) 

8.34 The problem of 8.3, except, = _ (^) ^ (y> t)) 2e^1 = _ ^ (y> t)> 

ap( lv 0,t) ~ A »oP (x, 0, t) = - (^) Vxo (x, t), and p (x, 6, t) = V»6 (y, t). The initial pressure 
p(x,y,0) =<fi(x,y) 

8.35 The problem of 8.3, except, - - (fc) *o (v. *), 2e ^ i +A V aP (a, V, «) - - (fc) V yo (y. fl. 

ap( ' lv 0 ' t> ~ A*oP(x, 0, j) = - Vio (x, «), and p (x, 6, t) = Vx6 (y, t). The initial pressure 
p(x,y,0) = v?(x,y) 

8.36 The problem of 8.3, except, = - (^) V» y0 (y, «), P (a, y, f) = W (yft), 

SE %T l£i " W(*,0,*) = -feW*,*), and ^E^l . -(#-Wx,t). The initial 
pressure v w/ 

p(x,y,0) = y>(x,y) 

8.37 The problem of 8.3, except, 5^ = - (^) V v0 (y, *), = - (fc) (y, fl, 

2E ^-A lO p(x,0,0 - -(£Wxo(M), and = -(fWx,0. The initial 

pressure v *' 

P(x,y,0) = <p(x,y). . . 

8.38 The problem of 8.3, except, <^ = - ^ (y, 0, 2E^2 +AyoP (a> „,*)--($) ^ (y. 0. 

- A l0P (x, 0, 0 = - (>) Vxo (x, 0, and fiEteM = - U) 0 l6 («, 0. The initial pres- 
sure \ »/ 

P(x,y,0) = <p(x,y) 

8.39 The problem of 8.3, except, ^E^il = _ ( jl) ^ yQ (y> f)) p (a> y> t) = ^ (y _ t)> 

~ A*oP (x, 0, () = - (,jL) ^ (x, t), and ^ + Aiftp {x> b,t)--U) ^ («, 0. 
The initial pressure ' 1 
P(x,y,0) = <p(x,y) ; 

8.40 The problem of 8.3, except, = - (f) ^ (y, 0, §E^1 = _ (^) ^ (y , 

- A l0 p (x, 0, t) = - (^) Vxo (x, t), and + Xxbp {x< b> t) ' m _ ^ ^ (x> t) 
The initial pressure p(x, y, 0) = <p(x, y) 

8.41 The problem of 8.3, except, = - (^) ^(y,t), 2E^!i+A vo p(a 1 y,0 = - (%) 1>y (»,*), 

- A x0 p (x, 0, 0 = - (£) ^» l0 (x, t), and 5^+A ltP (x, b, t) = - (#) ^ t (x,"i). The 
initial pressure p(x, y, 0) = v>(x, y) 
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~ Al0P (l ' 0> = " (#) ^ (*• *). and p (x, b, t) = Vx6 (y. t). The initial pressure 

p(x,y,0) = v>(x,y) 

843 

7^ " W (^0 ( 0 = -(^) <fco (x, t), and p (x, 6, t) = (y, t). The initial pressure 
p(x,y,0) = v>(x,y) 

8.44 The problem of 8.3, except, ^ - A^p(0,y,0 = - (^) ^'(y.i), 

"2* + A„«P (a, y t , — (fc) * vo (y, 0, &%2d _ W (1 , 0, 0 = - fe) ^ (x, *), and 
p (x, b, t) = rp xb (y, t). The initial pressure p(x, y, 0) = v?(x, y). . . . . 

8.45 The problem of 8.3, except, ^ _ AvOp(0>y>() = _ fej^t), P(a,y,t) - ^ (y ,' t) / 

ptSfe" AX ° P(X,M) = "(*)*•<»■'>• ^ = -(*)*.(..«)• The initial 
P(*,V,0) =v?(x,y) 

8.46 Theprob.emof8.3, except, -X v0P (0, y, t) = - (fcj^fcfcft^ - - (#■ ) to,*). 

s"^ ~ W (X> °' 0 = " (*) 0 *° ( *' * ^ - " (*) (*• 0 ^e iniJal pres- 
p(x,y,0) =<p(x,y) 

8.47 The problem of 8.3, except, _ Xy(>p (0> y>t) = _ (y _ ^ 

^ + A " p (a ' * f) - " (*) ^- (y, 0, - A I0P (x, 0, t ) - - ( t ) ^ (x , and 

P ' " ~ ^ x6 ( x> *)' The initial Pressure y, 0) = <p(x, y) 

*%T* - JUp(«,0, () = -(^(x,*), and ^ + A l6 p(x >M ) = -(f) **(«.«>. 
The initial pressure p(x, y, 0) = y>(x, y) Kv ' 

8.49 Theproblemof 8.3, except, ^l- XyQP ^ t) m _ '^i*tg*-'- U) M ,t), 
*V* ~ AxoP (x, 0, t) = - ( f.) ^ (., 0 , and ^ + Ax6p (x> M ) = - (f ) ^ (J, t) The 
initial pressure p(x,y,0) = <p(x,y) V vJ 

8.50 The problem of 8.3, except, - Ay0 p (0, y, <) = - (fcj ^ 0 "( y , ^ 

^ + A W (a, y, 0 = - (^) ^ (y , 0, ftJM - A^x, 0, t) = - (f) ^ (,,,), and 

+ <W (x, b, t) = - Vx6 (x, t). The initial pressure p(x, y, 0) = ^(x, y ) 

9 ^e a probl^" infinite lameUa rectangle ' incl «ding mixed boundary 

9.1 Cascadedsemi-infiniteIamella:atx = 0 1 p J (0, y ,0 = V' y i(y,t),6 i <y<6 i+1 , 

l^; 1 - ' > 0. At y = ^ = -(^) >( x, t , and at y = fc, 
5v H = ~ (^) N ^ H ( x > ')■ At the interface y = b h Vj = 1, K - 1 

<, (x, o = - = _ (te^j and 

A,^- (x,() = { Pi _! (x,6 i)t ) _ Pj . (x.ftj-.t)}. The initial pressure p^x.y.O) = ^(x,y) 
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9 2 Slo P „?) b,em °( 9 \ eXCept ' &t 1 = °' WC have a mbced boundary condition, which is 
hx = ~ Xfcj^vi (V.O for bj ,<y<b j+u j = 0, 1 jfc - l; 

Pi (0, V, t) = i> yj (y, t),bj<y< b j+1 , j = k , I- 1; 

= " (^) . (V,t) for 6, < y < = I, N — 1, {i>k+l}, t > 0. 

At y = bo, 2^ = _ (£) /Q ^ ftt y m ^ = _ At the 

interface y = Vi = 1, .... N - 1, ( I( *) = _ 0e) = ~ (^) • (^'ff' 6 '' 0 ) ' 

and AjV, (x, t) = (x, fc,, E ) - Pj . ( X) t )}. The initial pressure Pj . (x> y> 0 j"2 ^ y) 

93 l h 5o Pr S blem f 9 \' eXCePt ' 8t * = °* ™ have 8 mbced bou ndary condition, which is ' 
= " {t) j (V. 0 for bj < y < b j+lt j = 0, 1 jfe - 1, t > 0. 

Pi (0,y, 0 = Vvi (y,t), 6; < y < b j+1 , j = jfc, I- l, t > o. 

- A wPi (0, y, t) = - (y, t) for bj < y < b j+l , j = f, K - 1, t > 0. At 

y = &o, ^^f^ = - (£) o ^ 0 ( Iit ) and at y = 6k, P (x,&k,i) = ^ (x,t). At the interface 

J = Vj = 1 * - L ^(*.*) = - (fr)/*^) = - (*f) and 
\jrl>j (x,t) = { Pi _j (x,b jt t) - Pj (x,bj,t)}. The initial pressure Pj(i,y,0) = <pj(x,y) 

9.4 The problem of 9.1, except, at x = 0, we have a mixed boundary condition, which is 
Pi (0, y, t) = (y, t ) for bj < y < b j+1 , j = 0, 1, , it - 1, t > 0. 

- KiPi (0,y, (y, <). < y < j = k, r- 1, * > 0. 

2E ^ 2 = - (^)>w (y.O for *i < y < tyw, j = f, N - 1, t > 0. At y = 6o, 

P (x, 6 0 , *) « (x, <) and at y = 6 H) p (x, 6h, = ^ (x, 0- At the interface y = b j , 

= 1 N " * 0 - " ft), t 20 ^) = " (^#^) and 
AjVj (x, t) = {p^x (x, b jt t) - pj (x, bj, t)}. The initial pressure Pj (x, y, 0) = Vj (x, y) 

9.5 Cascaded rectangle: at x = 0, Pj (0, y, t) = Vw (y, 0 and at x = a, Pj (a, y, t ) I ^ . (y, t). 

• r 6j+1 ' ^ = 1 ^ t>0My = 6 °' 2E ^ = -fe) 0 ^ (*•') and at y = 6 N , 
^X"^ = - (#) r t/<K (x, «). At the interface y = 6,-, Vj = 1, N — 1, 

* M — (£) . (^i) = - fr)^ (^hJl) and 

Aj-^ (x, t) = { Pi _! (x, frj-, t) - pj (x, t)}. The initial pressure Pi (x, y, 0) = vj^^x, y) 

9.6 The problem of 9.5, except, we have a mixed boundary condition, which is- 
(») At x = 0, P j (0, y, t) = xl> yoj (y, t) and at x = a, Pi (o, y, t) = i> Vai (y, *), for 
bj<y<bj +1 ,j = 0,1 ,JE- l, t> o. 

' (") At x = 0, = - (^)^ ^. {Vt t) ^ at x = a, 2^ = - (fc) ^ (y> t) for 

bj<y< b j+1 , j = k, .: X-l,t>0. 3 

At y = b 0 , ^ = _ ^ (Xi t) and at y = 6h, ^ = - * H ( x , (,. At the 

interface y = 6,, = 1 K-l, ^ (x, t) = - (^) (^^'^) = - (V) ^-.jj.^.Oj 

and A^,- (x,«) = { Pi _, (x.dj.i) - Pi (x,6 JlE )}. The initial pressure Pj (x,y,0)= '^(x.y) 

10 Octant 
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10.1 The medium is bounded by the planes x = 0, y = 0 and z = 0 ; x, y and z extend to oo in the 
directions of x positive, y positive and z positive. At x = 0, p(0 y y t z,t) = 0, y > 0, z > 0 
t > 0, at y = 0, p(x,0,z,*) = 0, x > 0, z > 0, t > 0 and at z = 0, p (x, y, 0, t) f = 0,' x > o| 
y > 0, * > 0. p tend to zero as x oo, y oo and z oo. The initial pressure 
P(*,y,z,0) = 0 

10.2 The problem of 10.1, except for all t > 0, p (0, y, z, *) = (y, z, 0, y > 0, z > 0, p*(x, 0, z, tj = 
0x* (x, z, *) x > 0, z > 0, and p (x f y, 0, 0 = ip xy (x, y, t), x > 0, y > 0. The initial pressure 
P(x,y,z,0) = y>(x,y,z). y?(x,y,z) and its derivative tend to zero as x -> oo, y -* oo and 

£—►00 

10.3 The problem of 10.2, except the initial pressure p(x,y,z,0) = p/, a constant for all x > 0 
y > 0, and z > 0 

10.4 The problem of 10.1, except for all t > 0, p (0, y, z, t j = V y * (y, 0, V > 0, z > 0, p (x, 0, z, t) = 

(W), x > 0, z > 0, and M*fl0,0 = _ ^ ^ (xyt)j x > 0 , y > 0 . The initial 

pressure p (x, y, z, 0) = <p(x t y, z). y?(x, y, z) and its derivative tend to zero as x oo, y -+ 00 
and z — ♦ 00 

10.5 The problem of 10.4, except the initial pressure p(x,y,z,0) = p 7 , a constant for all x > 0 
y > 0, and z > 0 

10.6 The problem of 10. 1, except for all t > 0, p (0, y, z, t) = tf yx (y, z, y > 0, z > 0, p (x, 0, z, t) = ' 
^xx (*, *, «), x > 0, z > 0, and - A xv p(x, y, 0, t) = - (^) ^ y («, V, t), « > 0, y > 0. 
The initial pressure p(x,y,z,0) = <p(x t y t z). <p{z t y t z) and its derivative tend to zero as 
x — ► oo, y — ► 00 and z — ► 00 

10.7 The problem of 10.6, except the initial pressure p (x, y, z, 0) = p 7 , a constant for all x > 0, 
y > 0, and z > 0 

11 Infinite Circular Cylinder 

11.1 The medium is bounded by the cylinder r = r fl and extend to 00 inthe direction of r positive 
and (-00 < z < 00). The initial pressure p (r, 0, z, 0) = 0 

11.2 The medium is bounded by the cylinder r = r 0 and extend to 00 inthe direction of r positive 
and (0 < z < h). No-flow boundaries at z = 0 and z = h. The initial pressure p (r, 0, z, 0) = 0. 

A Special Functions 

A.l Gamma Function 

B A Table of Integrals 

C General properties and a table of Laplace transforms 
D Series 
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Octant 



10.1 The medium is bounded by the planes x = 0, y = 0 and z = 0 ; x, y and z extend to oo 
in the directions of x positive, y positive and z positive. At x = 0, p (0, y,z,t) = 0, y > 0, 
z > 0, t > 0, at y = 0, p(x,0,z,i) = 0, x > 0, z > 0, t > 0 and at z = 0, p(x,y,0,t) = 0, 
x > 0, y > 0, t > 0. p tend to zero as x — ► oo, y — ♦ oo and z — ► oo. The initial pressure 
p(x,y,z,0) = 0. 

A quantity Q of fluid suddenly injected at a point (x 0 ,yo>*o), [xq > 0, y 0 > 0, zq > 0], at time t = * 0 , 
[*o > 0]i and the resulting pressure disturbance left to diffuse through the homogeneous Octant. The medium 
is bounded by the planes x = 0, y = 0, z = 0,and extended to infinity in the direction of x positive, y positive 
and z positive. The pressure at the planes x = 0, y = 0 and z = 0 being kept at zero. The initial pressure 
is zero. 

We find p from the partial differential equation. 

dp d 2 p d 2 p d 2 v Q 

di = v *d£ + i y d£ + n *a? + ~ Xo)(J(y " yo)<J(z " ~ ' o) (1011) 

where rj x = t/„ = ^ and r; z = ^ with initial condition p(x,y,z,0) = 0, for x > 0, y > 0, 

z > 0 and boundary conditions p(0, y, z, *) = p(oo, y, z, t) = 0, p(x, 0, z, t) = p(x, oo, z, t) = 0, p(x, y, 0, t) = 
p(x,y, oo, t) = 0, t > 0. 

We apply the Laplace transformation to equation (6.1.1). We get, 
d 2 p d 2 p d 2 v Q 

t'dx* + ^dy 2 +r}x d?' sp=1 fa 5 ( x ~ Xo)5(y ~ yo)(J(z " Zo)e t 10 - 1 - 2 ) 

oo 

where p = / pe'^dt. We now apply the appropriate Fourier transformations to equation (10.1.2). We get, 
o 

1= gsin(nx 0 )sin(myo)sin(/zo) tn 

P <t>ct(n 2 ri x + m'ffc, + l 2 m + 5) llU ' 1 ^ 

where p = J ps\n(nx)dx i p - J psin(my)dy and p = /psin(iz)dz. 
00 0 

Successive inverse Fourier transforms of equation (10.1.3) yields 



P - 1 



f (*-x,)' ■ (y-yp)' ■ (z-xo)O * 
\ »7» T »?» T / 
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I *ht V* J \ V* + Vv + V. / 

\ »?« T »Jy I?, J 



p. 



Inverse Laplace transform of equation (10.1.4) yields 



(10.1.4) 



P = 



U(t- t Q )Q 



80ct{7r(t-t o )} 3 y/WhV. 
x < e 4 n*(*-*o) — e *o> > 



(10.1.5) 



where (* - * 0 ) = j J * > *° is the Heaviside's Unit step function. As t * 0 , pressure from equation 
(10.1.5) tends to zero at all points except at [x 0 , y 0 , zq] where it becomes infinite. 

With a change of notation the equivalent solutions for heat and mass diffusion problems may be formally 
written as * 



Heat diffusion: 



Qe- 3to s$ 
<t>ct (2irC p )% y/K x K y K z 



\ ** «V ** J 

l "» "» " « J I J 

(l£±52ii + lE±Soi2 + I £Z £oii\ i J(»-»o)» , , («+*o)Q * + 

I *• "» «■ J \ K« T <C» K, / 



'To convert from pressure diffusion to heat diffusion, we take t, = «, K = £,and ** = C p ; to convert to mass diffusion we 
take i} = £> and set tfc t = 1. * 
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Jlg+*o) 2 , (v+yo) a . (z±z^\ < 

\ K x ~ K y ^ K M J 



(10.1.6) 



and 



v = 



^wt-^wt e ^7\ e 

where the Q is the total quantity of heat instantaneously liberated at time r = t 0 at a point [x 0 ,y 0 , zo}. 
Mass diffusion: 



c = 



Qe- at °s$ 



(2tt)^ y/D x D y D z 



\ O, Dy D m J 



r/^« r~ it? r — -i — - + 



f (*+*o) 3 , (B-t>o) a , (*-*o) a 1 



Jl*±*o) a . (y+y 0 )» . [5+3^ \* 



(10.1.8) 



and 



c - 



«'{.-**-.-**} (I „,. 9) 

Here Q represents the amount of substance deposited at time t = t 0 at a point [x 0 , y 0 , zq). 
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The continuous constant source solution may be obtained by integrating the instantaneous line source so- 
lution with respect to time. However, for illustrative purposes, here we solve the problem in a formal way. 
Fluid is produced at the rate of q(t) per unit time from t = t 0 to t = t at the point [x 0 , y 0 , zq). We find p 
from the partial differential equation 

dp &p i d 2 p d 2 p „, x q(t-to) 

with initial condition p(x,y,z,0) = 0, for x > 0, y > 0, z > 0 and boundary conditions p(0,y,*,t) = 
p(oo,y,M) = 0, p(a?,0,z,t) = p(x,oo,z,t) = 0, p(x,i/,0,*) = p(x,y,oo,*) = 0, t > 0. Following the 
procedure outlined earlier, the successive application of Laplace and Fourier transforms reduces equation 
(10.1.10) to an algebraic equation, which is 

I _ Q (*) sin(nxo) sin(my Q ) sm(lzp) _ stQ 



<j>Ct(n 2 T} x + m 2 7} y + l 2 7] z + s) 



(10.1.11) 



Successive inverse Fourier and Laplace transforms yield 
q(s) e- st °si 



P = 



fat (2tt)* ^r} x 7] y r] z 



-X 0 )3 



(y-vo) 2 a (f- 



/ (*+xo)2 (y- w )a (z-* 0 )M * f (x-x 0 )* , (y+y 0 )* , (z-zp)M * + 

I J** Vv ^ V* ) \ »7x *7y »7* J 

J (x+* 0 )* (s- Zo )2 1 j f (x-x 0 )a . ( v -y 0 )* (z+z Q )a\7 

I 1« / \ »7« ^ riy ^ J 



f (x+x 0 )2 ( t ,- yo )2 f (x-xp)» 
I *J* fv 1« / \ »?» 



+ (y+yo) a + 



/(X+Xp)3 ( V +yo) a . (f+go}M< 
\ fas Ify Tig J 



(10.1.12) 



and 



p = 



- to) 



t-t 0 

/ gft-tp -r) f .(«-*q) 8 («+«n) a \ f _<»-»o) a (y+yo) a 1 
h J rl I J \ J 



8<t>ctn% ,/VxVvVz - 

f <«-'Q) a <*-Mn) a 1 

x < e - e S dr 

If g (t) is a constant and equal to 9, integration of equation (10.1.13) gives 

'r{*.^(^ + ^ + ^)} 



(10.1.13) 



P = 



»o) 3 ■ (,-»)» (,-«,)» 1 
1» T "fa J 
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r/l ; 1 f (*+x 0 ) a , (y+vo) 2 , ( £± zn^\ 1 



(*+*o) a , (y+vo) 2 ■ (z+z 0 ) 2 \ 1 
-v- -»/ y *7x T i» t?« y J 

V i ^ ^» t« / 



For the case where 9 (*) = g^, 1/ > 0, t > 0, integration of equation (10.1.13) gives 



(10.1.14) 



P = 



qr(v + l)U{t-to)t("-i) 
^4>Ct^ y/r\ x 7\ y r\ z 

VI 1- % •»» / 

// (*+*o) 3 1 (y-w)» ■ (*-*o) 3 ] 

V i % / 

l( {x+Xo)' (y+y 0 )* (z-Zp)M 

V I *»- 1y T 17, / 

/f (»-»o) 3 ■ (»-«>)» , (z+*n^7 + 

V l f« T «?y T »?. J 

lf (x+z 0 )> . (y-y 0 )» + 

y\ «?. T % T »», / 
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yp^xo)» | ( „+y 0 )» | (»y°) 3 j 



(10.1.15) 



The solution corresponding to the case where instantaneous and continuous multiple point sources at 
frot. Vol, zoi] at times t = fa, t = 1, 2, , JV, may be obtained by solving the partial differential equations 

(Pp d 2 p d 2 p 1 ^ 



dp _ d 2 p 8 2 p d 2 p 1 A 

dt~ r,x lti 2+ ^dy 2 + V'dz 2 + ^ t ^ Q ^ X ~ XOl ^ V ~ y ° l ^ Z ~ 201)5(4 " * 0t) (10.1.16) 



and 



dp _ d 2 p 8 2 p 8 2 p 1 A 

at ~ Vx d!? + ^dy 2 + ''a? + ^ 2^ y ( * " * 0t) Qi & ~ toi W x ~ Xot ^(f ~ f°') d ( 2 - **) (io-i.i7) 

respectively. The solutions for instantaneous multiple point sources in an Octant are 



V = 



\ V* Vy Vz J 



'1 GA^g" 2 


I (y- 1/0 *) 3 




((z+soi) 3 , 


(y-yoJ 3 , 
*?y T 


(*-*o*)M* 










(y+wu) 3 , 

>?y 


(*-*o.)M* 


(^/{^ 


1 (y-yoJ 3 

^y 






(y-yo.) 3 . 

*?y 


(*+*o) 3 1 * 


*l (\/{^ 


, (y+yo*) 3 

»7y 





( (g-soi) 3 , (y+yo t ) 3 , (z-r 0t ) 3 1 3 

\ T/x T tj y V* ) 



{ (X-XQ, 



*7y 



*?y 



T7* 



and 



P = 



1 »7* *?y *?a J 

1 f HLJoja f- ^ 



(10.1.18) 



f -JilzWjJL (y+ytu) 3 I 

^-'oTT > < e ^v^-'oJ — e 4f iy(«-*o*) J. 



(10.1.19) 



1 
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The solutions for continuous multiple point sources are 



P = 



si 


N 










f(x-x 0 i) a . (ir-yo 4 )a (z-z<u) a 1 * 




»7y 1 »7« J 





I * ^ J I »7* + *?y + il. J 

f (x+x,,)' (g+jffiti; + i*^ig i il\* f (x-x 0t V» (y-itt.)* , («+«».)» 1* + 

I * ^» •>» J In. T ♦»» "•" t(. / 

I *>■ "» * J I «?* + r, v + tj. / 



and 



f (x+soQ3 (y+y 0t ) a ■ (s+zoQ' V 
I T »?y T 1. J 



t-to* 



(10.1.20) 



p = - 

8(/>C t TT 



te^j ( 0t) / — n — r ^^- e -^}{e -^- e h^| x 

e S#-_ e -^-j dr (1Q121) 

We consider some special cases of practical relevance. 

(i) A line of finite length [zq 2 - z 0 \] passing through (x 0 ,y 0 ). 

The solution is obtained by simple integration. For an instantaneous source, we get 



and 



U(t-t 0 )Q ( -.<,*-*q)\ _J*+*a2Ll f (v-^o) 2 ^ <g+m>i 1 
S^K^^^ r WHT - e *^}{e WT-e J 
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* {«" (iro) + ^ (to,) - - - " f (535^5) } 

(10.1.23) 

Solution for a continuous source is* 



and 



"W^Mw)--(l^)--(33)}* 

(10.1.25) 

The spatial average pressure response of the line \z 02 - 201] is obtained by a further integration t. 



+ 

and 

t-t 0 



-y — I e "* T +e i.r j | e ^T^-e ^ e _ e ^n^l dr (10.1.27) 



*The solution is relavent to pressure transient methods associated with partially penetrating or horizontal wells in the fields 
of ground water movement and the production of oil and gas from hydrocarbon reservoirs. 

txhe general results given by the equations (10.1.24) and (10.1.25) assumes that the pressure flux is uniform along the well. 
A good approximation to the wellbore pressure is obtaind by computing the spatial average pressure along the line. 
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(«) Multiple lines of finite lengths [zo 2l - zon], [x 02 t - *ou] and [y Q2 , - y 0 n) passing through (x 0 i, yot) for 
l = 1,2...,L, (y 0t ,2ot) for t = Z, + 1,2...,M, and (x 04 ,*0t) for i = M + 1,2..., N respectively. Where 
{L<M<N)\ 



P = 



L 



1 ^ 



and 



"W^) + -(^)X^)--(w)}* + 



*The solution corresponds to the case where there are sets of partially penetrating verticlal wells and horizontal wells in an 
artesian aquifer or hydrocarbon reservoir. 
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*MwW(wH(wH(w)W 

The spatial average pressure response of the line [2020 - *oioL * = 0, is obtained by a further integration. 
1 



P = 



$3 (5) e-* t0 > J ~p {cos(ti2oio) - cos(u2 0 2o)} {cos(uzon) - cos(uz 0 2t)} x 

+ ! x 

<£Ct7T 2 (^020 - ^oio) yfWh 

x_E ft M e-oJ?!£M {c os(uxon) - cos^)}/^ {J{^f^ + + " 

. 1 

0Ct7r 2 (2020 - *0io) y/Wh 
x £ * (.) {cos(uyou) _ cos^}/^ {^{(£^ + <£^}(. + , S )J _ 



and 

P = 



80ct7r (2020 - 2010) v'^M/v 
x * 



v^ rm f x / Qi(t-t-u) f _i£-oji ,<"j-Q>>M f -<yoi> a .(v+Ulklil 
/ „ (t — tot ; y ^ e — e J | 6 " — e 4,?v | x 
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+ 
+ 



+ (^ + «..)«f(^^)-(^ + ,..)«f(2^) 

+2 yM _ e - i- c ...)' j + 2 yM|.- '-"O u) ' _ e- l- "Cr-H . 

- (% , + am)erf (2^)_ (wt ^ )erf (£^) + 

+20 { e -MS^ . e - } +20 {«- ^-O^ - J] * + 

t — tot 

* {- ( w) - " f ( w) - - ( w) + ^ ( w) } " 

4. I lilt + \ f ^(t-tOt-u) f -<*.-*P>)' (*+«n.) 3 1 

V^-w^jfe, ( } / — vs — r r 
x ( erf ( w) - f (w) " - ( w) + - ( W)} * 

* Hte) -<m) 

The solution of equation (10.1.2) in the infinite region [-00 < x < oo], [-00 < y < oo] and [-oo < z < oo] 
may be obtained by taking the complex Fourier transform [equation (2.2.1)]. 

= Q e inxo e imyo e ilzo e -8t 0 

P = <frt(n 2 T, x +m>T, y +l*T, z +s) (10 - 132 ) 
Successive inverse Fourier* and Laplace transforms of equation (10.1.14) are 



(10.1.33) 



and 



! . (x -.* n)2 . I ■ («-'n> 3 1 

U(t-to)Qe \ 4 "*<'-'<>> + <■>»<'-«<>> + <-7«c-«o> J 
P = Tn~, n4 (10.1.34) 



*See equation (2.2.11). 
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respectively*. The corresponding continuous solution is given by 



and 



T _ U(t- to) 7° g (t-tp- r) -{^+*&+*Hg£} 



We consider some special cases of practical relevance, 
(i) g (£) is a constant and equal to g 



and 



V = 



(«) g(<) = 9«", ^ >0, <>0 



-*o) a , (y-yo) a . (*-*o)M 

»/* *?v »?* / 



and 



P = 



_ _ gr („ + 1) ^ { tod! + tori! + fcal! } 5 ) 

/ i\ qj ' i ( < x - x n) a i («-»n> 2 , <«-'n) a 1 



(10.1.35) 



(10.1.36) 



(10.1.37) 



(10.1.38) 



(10.1.39) 



(m) A line of finite length [zq 2 - z 01 ] passing through (x 0 ,2/o). Instantaneous source 



2to?r 2 v^%. 



(10.1.41) 

•For translation into cylindrical coordinates we write x = rcos0,y = rsin0,x o = rocosflo.yo = r o sin0 o and tj x = rj u = rtr 
Then for an instantaneous point sourc equation (10.1.16) becomes 

f {r a + ra-2rr 0 co«(g-g 0 )> , z _, n >2 ] 

^— / ^^e L 



and for a continuous point source equation (10.1.18) becomes 
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and 



P = 



i f («-» ft )2 (y-vn) 2 1 

&4ct*{(t-to)}^ \ \2y/v.(t-to)J \2y/vAt-to))l 1 ' 



and for a continuous source 
g(s)e- ato 



P = 



{Ji^^WY^} [sin {/ {z ~ - sin { ' {z - *° 2)}1 dl 

(10.1.43) 

and 



For the case q (t) is a constant and equal to q 



P = 



qe -st 0 



u 

(10.1.45) 

and 

, . ^=^?i{ - (^. - f«=ai)\.-*{'* £ ^P e L ^ 

(iv) Spatial average pressure response of the line [202 - 201] 

(10.1.47) 



and 



I J 1 " V Zv/^tZ / \ 1TU\ ZQ2- ZQ! ji 

(10.1.48) 

(u) Multiple lines of finite lengths [zo 2t - zon]. [z<)2 t - xon] and [y 02l - yon] passing through (x 0l , y 0l ) for 
1 = 1,2...,L, (ycu.-SOf.) for t = L + 1,2...,M, and (x 0t ,zo t ) for 4 = M + 1,2..., N respectively. Where 
(L < M < N). 



P 24>c t 
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x [sin {u (z - *on)} - sin {u (z - zq2 L )}] du + 



+ —4— £ fc (.).-* 7iif 0 /^{(izi!!j! + k^\ ( . +tt ^)\ x 

x [sin {u (x — xon)} - sin {u (x — X021)}] du + 

x [sin {u (y - y QU )} - sin {u (y - y 0 2i)}] d« (10.1.49) 

and 

1 ^ 

S(f>ctn y /7} x 'qy 

l 

_j — z= x 



X*~^ 
=L+1 0 



^-«T tfi= ^{-(^)--(^)}-'* { "* e+fc * fi} ** 



0 

(10.1.50) 



: t r(( ^;/^^ {etf( ^). erf(w)}e -*{^^L 



The spatial average pressure response of the line [zo20 - zoio]> 1 = 0, is obtained by a further integration. 
tyctiPyfrg^ ^ (2020 - *oio) 

cx> 

x J u \l cos i u ( Z <>10 ~ ^1*)} ~ cos ( w (*02$ - son)} - cos {u (2010 - *02t)} + cos {u (Z 0 20 - S02t)}] x 

0 

Z<t>W 2 y/rh% (zo20 - S010) 

OO *02O ✓ / x 

: yi [sin {u (x - xon)} - sin {» (x - x 02t )}] J K o\J{ + } (» + u* Vx ) \dzdu + 

0 ZniA * * ' 



0 roio 
2^>C t 7T 2 N /^^ t J^j (Z()20 - ^010 ) 
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x J - [sin {« (y - y ou )} - sin {u (» - y 02l )}} J K 0 j ^{^^ + } (* + u%) \dzdu 



(10.1.51) 

and 



V = 



8</>Ctn (2020 - ^10) y/n^h 
t-t 0l 



.0 /^*/ -^°n" a ° u > a ('Qio-'ou) 2 ('o^-'oaJ 2 ifoiozfo^l -Al ^^^ + ^'^^ ll 

V~t e ~ e ~ e +e r e * J du+ 

*MwH(w)} e ^» <—> 

10.2 The problem of 10.1, except for all t > 0, p(0,y,z,t) = if> ys (y, z,t), y > 0, z > 0, 
p(x,0,z,t) = t/>* x (x,z,t), x > 0, z > 0, and p(x,y,0,t) = i/>*y (x,y,t), x > 0, y > 0. 
The initial pressure p (x, y, z, 0) = v?(x,y,z). <p(x,y,z) and its derivative tend to zero as 
x — * 00, y — ► 00 and z — ► 00. 

The solution for the continuous point source is 

I T Vy T »7. J 

/l£±xol! , (y- W )» , (z-zp)» 1 i f (»-x„)» („+„„)» (z-z„V» ) * + 

J (x+» 0 )' (z^zp)' ]* f (*-» 0 )' i (y-yo)' , (z+z„VI* + 

I *** 1» 1« J \ >»* + «7» + IJ. J 
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+5 



f (*+*o) a , (V+Vo) 7 , (*+*o) 2 1 * 
\ *fc ^ »?v "■"»?, / 

nr )xi> yz {m, I, s) + mr\$ xz ( 




ooo 



n, /, s)+lr} z ip xy (n,m, s) l 
,m 2 + *M 2 ) j 



^ woo oo 



sin(nx) sin(my) sin^dndrmtf + 



; -uQ 2 j ' 



f (*+tx) 2 . (y-v) 2 ■ (z-w)*\< 



i » * _ i« / \ -?» + %^ + Nr~j 

I *» / I >7x + r,„ + 



f (*+u) 2 . (y+t,) 2 , (z+tx;) 2 !* 
\ n, T t? v * T ^ J 



dudvdw 



(10.2.1) 



and 



V = 



U(t-to) 



t-t 0 



/ V(t-to - T) ( _ (x-x Q )3 (x|» 0 )3 1 f (i,-vn) 2 (y+y n ) 2 ^ 
1 < e 4t »« t - e We 4f >v T - e 4l »w T i 

v - - Q T * I / I / 

f («-*n) 2 (*+*n) 2 ) 

^ oo oo oo t 

+ A / / / / { nr >^>vz K /, * - t) + mJ7„^ z (n, Z,< - t) + Zt; J iy (n, m, < - r) } x 

n n r\ e\ * 



0 0 0 0 

x sin(nx)e- ,J * naT sin{my)e-"y m ' T sm{lz)e-*' l * T dTdndtndl + 

OO OO OO 



8 
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t 

I 

I 

I 



p 

and 



x{e-^- e -^} dudvdw (io22) 

y(n,m,s) = JJi> xy ( X ,y, S )sMruc)sMrny) dxd y, ?.,(„,,,,) = / jL (x, 4 ) sin (nx) sin (/z) 
^ KM) = J{M'**)*l«»)*>M+iM, ^(„ >mii , = / jl s (x >y ,Osin(nx)sin( mi ,) da:dj/) 

We consider some special cases. 

(0 (y, z, *) = *y, (t) , (x, *, f ) = ^ («) , ^ (*, y, t) = ^ W) 9 (f) = 0 md ^ V) w) = 0 

= ± 777f n2 ^vz (s) + m 2 V»* xz (s) + p v rf n ( S ) ) 

* 3 J J J \ nml(s + Ti x n* + Vy m 2 + Vz p) j ■"»(«*) sin(my) sin(Zz)dnd77wfl (10.2.3) 

^_y_ f ip xz (t-r)e *<»r .( X \ / z \ 

+ ^ / ^ (< "v5 C - CTf (2^7) erf (2^) * J dO.2.4) 
^v,w)=0 )=PVt ' ^ (x ' M) = P "' <«• *> *> = P*v-' P** Pxv are constants. 9 «) = 0and 

OO OO OO 

n3 { { { \ snml (* + + f) v m 2 + 7fe/2) / sin ( nx ) sm(my) sm(lz) d n dmd l (10.2.5) 

and' 

if 5- - ** 
t ^ 

(m) y, ,) = M. ; x > o, „ > 0, , > 0, g (0 = 0 and %z (y, *, ,) = ^„ z , t) = ^ (l> y _ t) = Q 

f ^ r .((( sin(nx)sin(my)sin(fz) 

Pl JJJ(n*v* + m\+l2r, z + s) dndmdl (10-2.7) 
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and 

(»v) «,(«, y, z) = -^=; z > 0, y > 0, z > 0, q (t) = 0 and ^ (y, z, t) = ^„ («, t) = ^ (x, y, t) = 0. 
p = Pl J( 1 V 7 7 7 Sin(nl) SiD(my) Sin(fz) r r rr 

and 



p = 21(1)* f x 2 \ f y 2 \. / z 2 \ 

(v) A line of finite length [zq 2 - *oi] passing through (z 0 , y 0 ). 



f y m 2 + 77 2 / 2 ) ?sin(nx)sin(my)sin(/z)dndm«« + 




^ 00 00 00 

+ / 9 x§ 3 / : / / / v(«.w,t») 



(s + 77 x n 2 + T]yiri 



-«)» ■ l»-v)» , 



1 1* T 1y T »Ix / 



\ »7* T «7. / 



dudvdw 



(10.2.11) 
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and 

t-t 0 



v = U(t-tp) f° q(t-t 0 -T) ( i*=p£ <-+->» 1 f <.-»,»» (>+Bn) " > 

^oooooot v 7 ' 

"*W / / / i nV '^ (m ' 4 - T ) + (n, /, * - r) + Ir, J ly (n, m, < - r) } x 

0 0 0 0 J 

x stain*)*-*-"- siii( m „)e->."'' smlkje-'-'"'**^!! + 

1 i (10.2.12) 

The spatial average pressure response of the line [ ZQ2 - ^ „ obtained by a ^ 

x sin(nz) sin(my) {cos^j ) - cos(/z 02 )} dndmdl + 

+ tf / // (* + ") 2 (y + y)2 1 )1 

°1V1 + ~^/ (5 + /2% )]J sin (H{cos(Z^ 1 )-cos(^o2)}d/d u d V d U ; (10.2.13) 



and 

t-t 0 



: ^ ^ ^= j iiL^_l) J (2fl2 _ ^ erf (M-J^ + {ZQ2 + ^ erf 'f^, ^ + 
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^ oo oo oo t 

+ 7r 3( 202 _ 201 ) J J J J] { nr >J>v* ( m - M - r) + (n, I, t - t) + ZtjJ^ (n, m, * - t) } x 

0 0 0 0 

x sin(Tuc)e-"* n ' T sin(mif)e-"v m ' r {cos(iz 0 i) - cos(/zo 2 )} e-i'^drdndrndl + 

OO OO OO 

"{-(^) + -(^)-(i3ff)-(i3i3f)}*— 

(vi) Multiple lines of finite lengths [202, - zo U ], [x 02l - x ou ) and [y 02t - yon] passing through (x 0t ,y 0l ) for 
1 - l,2...,L, (y 0l ,2bt) for 1 = Z, + 1,2..., AT, and (xot.^oj for 1 = M + 1.2....JV respectively. Where 
(L < M < N). 

1 

P ~ ' X 

+- = X 

jv/ 00 r t i~ ■ — 



+ 1 
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oooooo ( ^ ^ 

*S £ I (s + »fcn2+ ^2+^/2) Wnx) sin(my) sin(iz)dndmd/ + 

^ oo oo oo 

— -|-— ZZ / / / <P(u,v,w) 



dudvdw 



(10.2.15) 



and 



+ £7 ■ V) - fp.) / ft (* - «fc - «) / -ferfg - (^n.)' \f 

S^vW^i 7 u \ " e je _ e ~5iit- I x 

^ oo oo oo t v 

(m ' * - T ) + («. M - r) + ft,,? („, m , t _ T) l x 

0 0 0 0 'J 

x sinCnxJe-"^ 3 - sm(my)e~^ sm(lz)e-"^ drdndmdl + 
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oo oo oo 

%(^!^///^ (u,u,w) { e "^- e "^}{ e "^- e "^} x 

e ^-e-^-jdudvdu; (1Q216) 
The spatial average pressure response of the line - ^j, 4 = <>, is obtained by a further integration. 



x X> (*) e -,t0 * ^ ^ {cos(«2oi«) - cos(uzo2<>)} {cos(uzon) - cos(u2 0 2i)} x 

1-1 0 

+ 1 JJ 

-*{i/{^^}c^>} - * {^^^F^} + 

+ 1 

fa* 2 (2020 ~ ^oio) y/Wh 

000000/= = = . 

I 8 f f [-1 nT k^»* ( m ' *' s ) + TO %^xz (n, f , s) + Z^ xu (n, m, a) 1 

t 3 (2020 - A)io) J J J I I (* + »J.n 2 +»J B m a +i| I P) f * 

x sin(nx) sin(my) {cos(/2 0 io) - cos(izo 2 o)} dndmdl + 
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-K 0 
+ K 0 



° | ]j{ + } Z + < 2 ^) j j nn(to) {cos(/zoio) - oos(^ 2<> )} e^'dli 



l dldudvdw 

(10.2.17) 



and 



1 

— x 
8^*" (2020 ~ *Ol<>) \/^7y 

+( »» + w-(a^)-h» + *)-(^i) + 

- <~ - -» - (^r ) - <- - - 1 52 ^ ) - 

- h - + * J -'(^ t )-(*» + »J-( a ^ t ) + 

1 Af t ~ 1 o* 

- {- ( w) " - (^?) " - (^0) + - (^) } « 
"{-(^f)-(^»)-(^)-(^)}*. + 

1 AT *-*0» 

^(^o-zoio)^^! / ^ 1 ^^- e ^^} : 

^ OO OO OO t 

0 0 0 0 



'*7z^xy (n,m,t-r)|x 
cos(/zq2o)} e" r, ' l2T dTdndmdl + 
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^ oo oo oo 

+ «"(*0-*.0)VK*// / {'"^ - '"^J {e-^ -e-^l x 

10.3 The problem of 10.2, except the initial pressure p fx t/ z 0^ - « * «. ■ 

y > 0, and * > 0. p*«*»ure p ^x, y, z, Oj - p 7 , a constant for all x > 0, 

The solution for the continuous point source is 

5 



**(\/{^+ t 'y ; + i ^}*) 



+ 4 

7T 




ggx^yz (m, /, 5) 4- mr) y i) X2 ( 

(5 + VxTl 2 + fy™ 



w »^ 

— f f f sin (' u; )sin(m3/)sin(/z) 
* 3 / / / l« + V*n* + T, y tn2 + „,/») nml dndmdl 



[n,l,s) + lr lt ip sy (n,m,s) ) 
t m 2 + Vx P) j 



0 0 0 



sin(m) sin(my) sin(/z)<fndm<i/ + 

(10.3.1) 



and 



P = 



g oo oo oo t 

+ ^ II II (m ' « - T ) + ™J^« (». * - r) + 6|,^_ (n, m, t - r)l 

0 0 0 0 J 



62 



— oo oo 

^^K^^Jp..(-.^)^(»)^ (l2) ^j„ K/ ,, )= 77^ (ti „ )8to(m9)sto((z) ^ i 
*..<«.»■..> = 7h^,,,^^ (m , )dx „, 5„ ( „, M) . 77L (I ,», !)8in(ra)slI1(j2) ^ 

— OO OO 0 0 

(..«>- h.. <.,*, o („ ., . ]>„ „ e -., tU °2 i ^ <*,„,.) =7, Ml , i ,, ()e -.. (i , 

We consider some special cases. ° 
(») A line of finite length [zo 2 - zoi] passing through (x 0 ,y 0 ). 



f SiM {cos(l2Ql) _ cos(/ , 02)} |^|&za)! + (s + ~| _ 

,m 2 + ^/2) > sin (^)sin(m2/)sin^ 




™7s^ (ro, /, a ) + ro??^ (■ 



VJU OO 

f f f si n(^)sin(my)sin(^) 
^ / / / I* + + ^ + S?) nml dndmdl 



(10.3.3) 



and 



t— to 

^ oo oo oo t v ' ' 

^ m *~ K l,t ~ T)+ mr> ^" (ra ' z ' * ~ r) + ^ (n ' m - 4 - T) } x 

x sin(nx)e-"^ 3 - sin(my) e -^ sxn(lz) e -^ drdndmdl + 

^(W)^*)-^) (10-3.4) 
The spatial average pressure response of the line [zo2 - ^ is obtained by a further integration. 

p = (Z02 - a*) ^/pW^) - cos(/, 02) } 2 jjroj^l&zaJ! + kzae J (s + ^ J _ 
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+^rr-^ J 77 -I nr) *$y* (m ' f ' g ) + rggJW, (n, /, s) + fyj (n, m, s) 1 

x sin(nx) sin(my) {cos(fcoi) - cos(/zo2)} dndmdl + 

+—*EL 7 7 / ^Msin^fcosq^) -cosfi^ } 

(*02 - zoi) y 7 J {s + 7j x n 2 + T) y m* + r, z P)nml* dndmdl (10-3.5) 



and 

t-fo 

P = 



4^ (2 - zol) ^ / £iL V~ jl [<«» " ^)erf (^^) + (z 02 + *oi)erf (^&) + 

oo oo oo t 

^(AB-ibi) II f/l { nr,x *v* {m > * - r ) + ™&$» (». '> * - r) + fti,^ (n, m, t - r)\ x 
oooo J 

x sin(nx)e-'-^ sin^e-"™*' _ e -^r dTdndmdl + 

(10.3.6) 

f - ^ °/ T [ To" Z ° U] ' ( *° 2 < - * 01 ' ] and ^ - »» J ^ trough (z 0t , y0t ) for 



1 

P = — — . x 
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g) + (n, J, 5) + ( n> m, s) ) . . 
(* + ihn' + vrf+%P) jsin(nx) sm(my) sm(lz)dndmdl + 

, 8p/ 7 / / _sin(nx) sin(my) sin(/2;) , J Jf 
^ / / / (5 + ^ + ^ m 2 +r?z/2)nm/ ^ M (10.3.7) 




+ JL / /Y J nT k^»* ( m > s ) + m Vy^xz ( 



0 0 0 

00 00 00 



_ 1 / ^ f li{t~ t 0l . - tt) f _<*- * Bt ) 3 _ (*+*o t ) 2 "I f -(«-»0,) 2 . (v+vn.) a 1 

^*y/wh> 'J « t ^^" e ^^}{ e ^^-e "^-jx 

''{-'(^») + -(^f)-(^f)-(^)}*H- 

+ I stttu qg* - fru -«) f -ii^ai .(.j-^n r _<»-yfi.> a _(»+»fl,) 3 i 

"WwM^)--(^)--(i^)W 

^'^jfe-! ' 0t) 7 » \ e ^^-e 1 ^||e -c "#jx 

" Mw) -(w) - ( w) W 

00 00 00 t 

n, m, i — r) > x 

0000 

x sin(m)e- ,J * nJT sin(mj,)e- , '>' maT sin^e-"''' 7 drdndmdl + 
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The spatial average pressure response of the line [z* 2<> - z QU> l t = 0, is obtained by a further integration. 
1 



fan 2 (*02O - 2010 ) y/Wht 



7 l 

( 5 ) e " st0t / ~2 {«*(*«bio) - cos(u2 02 o)} {cos(uzon) - cos(uz 02 J} x 

H 2 _____ x 

^ct?r 2 (zq20 - 2^10) ^/t^ 

j> W e-oJ^fE) {COS(UIOU) _ cos( ^ 2t) ;jL 0 |^(i^ + (^j (s + u27?x) | _ 

{ \/{^^}<— )} - * { /{SSE}^} + 
-WHf 2 ^}'— )}] — 



X 



^TT 2 (^20 ~ *01<>) 



. 8 f f fl) n ^yz (™>> l > s) + rnr) y ^ xz (n, /, 5) + frfcVL, ( 

* 3 (2020-2010)/ J J I I (* + i7»n 2 +f7 y m 2 + ^P) f X 

x sin(nx) sin(my) {cos(/2 0 io) - cos(/2 0 2o)} dndmdl + 

+ Spi f f [ Mnx)sin(my) {cos(lz 010 ) - cos f/^U . J Jf 

7T 3 (^020 - zbio) J J J (s + Vxn*+ Vy m*+ Vz P)nmP dndmdl t 10 ' 3 ' 9 ) 



000 

and 



P = 



8<to?r (zq20 - ^do) y/fgrfy 
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_ / 77,11 ( (*02A-*OlJ 2 ^ \3 "\ / , 



-Kw) 

+ Vv tB ^ 12 - «" tffl « ui! } + {. - - .-ana* } - 

- <~ - ~> erf te*) - <- - *•> - (W) - 

, AT *-jo» 

+3T-7 1 c V U(t-U\ f 9* (* ~ *0t ~ ") f -i2f£f >s 



*iert( Z TZ? 1 ) - erf f * 01 <> ~ **\ . er f f 3»»±M . , /*oiO + 1 
! ' / ^^)+erff^^^_erf| < £z£o2 i \_ CTf /x + xo 2l \1 

*-S7~ 7 1 U(t-tn) f ^ft"fat-«) f -<*r*iu) 2 
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10.4 The problem of 10.1, except for all t > 0. v(Q u z t\ - .h < 

The initial pressure p (x, y, z, 0) = v ( x vz) u>fc « \ IP t~ v . (X> V ' ^ X > °» » > 
x - ee, y - oo and z-> oo v(»,y,z) and ,ts derivative tend to zero 



The solution for the continuous point source L 




68 



and 

t-t 0 



P = 



f ('-«n) 2 («+«n) 3 1 

xje -fnfr- + e — 1#- l dT + 

oo oo oo t 

R 

+ 



HI 1 1 { nVx ^ (m ' ' - r ) + m? ^« M - r) + -^ xy (n, m, * - r) 1 x 
0 0 0 0 v* 0 * J 

x sin(m)e-""» 3T sin(my)e cos(/z) e -"«' , ^ r rfndm(« + 

x je- 1 ^- + e - *^" j dududto ^ 4 ^ 



where (n, i, ,) = J (*, s) sin (n*) cos (Iz) dxdz, (m, /, 5 ) = J J ( y , s) sin (my) cos (lz) dydz , 

= oo oo_ 0 0 

^ y (n,m,s) = fSrl> xy (x,y,s)sin(nx)sin(my)dxdy, * xz (n,l,t) = SfrP xz (x,z,t)sin(nx)cos(lz)dxdz, 

= oo oo 0 0 

i> yz (m,l,t) = jfi> yz (y,z,t)sm(my)cos(lz)dydz, ^ xy {n,m,t) = JJj xy (x>y>4) sin (nx) sin (my) dxdy, 
We consider some special cases. 

(0 * v , (y, *, 0 = Vyz (*) , tf« (», z, t) = 0« (t) , ^ ( x , y , t) = * xy (t), q ( t ) = o and „, «,) = 0. 



OO OO 

-*// 



I^jMfL , rnnJ*As) i> xy {s) e -^V^ ^ 



0 0 

x sin(nx) sm(my)dndm (10 4 3) 

and 

P = ~ [ ~ { ( * " r > erf f-Q 4. ft - r) ejg / x \ 

^Jvt\ r^- x erf iw; + — ^ — OTf (vp) + 

+ ^ CTf (2^) CTf (2^) K (10.4.4) 
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vlw)=0 )=PuZ ' ^ (X,2 '* )= P"- ft., ft. and 9xtf are constants. g(t) = 0and 

P = 4 / /j 7 n77 f + . 

*r 7 7 |m(. + ^ + ^) + n(s + Vx n> +Vy m>) 

x sin(nx) sin(my)dradm 

(10.4.5) 

and 

+ ^ (2^?) CTf (a^) ] dr (10.4.6) 

(m) *>(*,„,*) = ^Lj; x > o, „ > 0) , > 0, ,(«) = 0 and („,*,,) = ^ (x , 2ji) = ^ (a . ^ = Q 

p = p J(l) 3 777 sin(nx) S m(my)cos(/z) 

V W / / / Vnm/(n^ I + m2 , ! , + ^ + s) d ^ (10.4.7) 



and 



8 U; V (a*) 'i (s^J (g^j (10.4.8) 

(if) A line of finite length [z 02 - z 01 ] passing through (i 0l y 0 ). 



'/y J j 1 

000000/ = = 
1/ /7j n7 lx^yz K I, s) + mr] y ^ xz (n, /, g) + ^ 



+ -3 

bob 



^ 00 00 00 

+ , 9 "1 . ~ / / m«,w,t«) 



(s + r) x n 2 +Tj y m 2 +r) z P) >sin(nx) sin(my) cos(lz)dndmdl + 

Mm****?*'*?}-) 
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\ t?x *7y »?* J 



dudvdw 



(10.4.9) 



and 



t-t 0 

_ ^(* - *o) / q (t - t 0 - T) / _<*-*q> 3 c+.^a -| f (v- yn ) a (v+i/n) 2 ^ 

P " 80W^ / ; l e ^-e ^jx 

oo oo oo t 

J 1 1 { ni)x ^v* (m ' Z - f - T ) + m vL i, * - t) + -jL^ ( n , m, t - r) J x 
x sm(nx)e-** n3r sm(my)e-*" mlT <ns(lz)e-*' l ' T dTdndmdl + 

OO OO oo 

xje ^?it + e (10.4.10) 
The spatial average pressure response of the line [z Q2 - z Q1 ] is obtained by a further integration. 

I 8 7 7 7 1 f n V*%z (m, l, s) + mr, y ^ xz (n, I, s) + -0 xy (n, m, s) ) 

t 3 (2o2-2oi)7 J J l\ ' (s + Vxn 2 +r] y m 2 + ri z l 2 ) f* 

x sin(na:) sin(my) {sin(/zo 2 ) - sin(Zz 0 i)} dndmdl + 

OO OO OO OO r- ^ i 
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+ k J // ( x ± u ) 2 (y + v) 2 1 1 1 

IV \~nT~ + ~^T~ J (5 + PVz) j \ COs(/u; ) ( sin (^2) - sin(/2oi)} dWudt/Ai; (10.4.11) 

and 

rr/ \ 4~to 

g oo oo oo t 

+ * 3 (zo2-zox) JJJJI \ nT, ^y (m ' Z > ' - T > + (». I. * - r) + -^ x „ (n, m, * - r) 1 x 

0000 0Ct " 'J 

x sin(nx)e-"^ 3 - s in(m^"^ {sin(/zo 2 ) - sin(Z^)} e^drdndnuU + 
(v) A line of finite length [x 02 - x 01 ] passing through (y 0 , z 0 ). 



+± f f f) nJ]x%S) ** (m ' *' s ) + m7 ?v^x» (n, /, a) + (n, m, s ) } 

| (s + %n 2 + 77 v m 2 +»fe/2) ~ >sin(nx) sin(roy) cos{lz)dndmdl + 



+ (2* 



| (»-t») 3 + (y-v) 3 + 
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\ I* 1?B 1|* J 

and 

t-t 0 



dudvdti, (10.4.13) 



"M^M^)--(3j©)--(3©)}* + 

^ oo oo oo £ 

+ ^JJJ J { nVx ^ (m ' l ' t ~ T )+ m Vv^*z (n, l,t-r) + -L? („, m> « _ T ) } x 

0 0 0 0 ^ J 

x sin(nx)e-"»" aT sin^e""*"* 2 ' co8(k)e-»- l>T drd»ulmdl + 

x {e-^+e-^Jdud^ (104 14) 

The spatial average pressure response of the line [x 02 - * 01 ] is obtained by a further integration. 

OO 

+ 8 ///I l 7 " 7 *^* ( m - '> 5 ) + m ^xz (n, /, a) + xy (n, m, s )| 

^(ate-aroOy J J n| (• + U,n» + vrf + „,/») j * 

x cos(k) sin(roy) {cos(nx 0 i) - cos(nx 02 )} dndmdl + 
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" K * I ^{ ^ + ^ ^ } ( g + n2 ^) I j sin(nu) {cos(nx 0 i) - cos(nx 02 )} dndiKfodu; (10.4.15) 



and 

i-to 

p = ^(*-*o) / 9 (* - *o - r) f, \„c( x m~ x o\ 



w^^-gs?)-')--- (^)— m- 

^ oo oo oo t 

+ 7t3(x 02 8 -x 01 ) / / / / 1 {m ' l < 1 ~ T ) + Mt-r) + j^ xy (n, m, * - r) } x 

0 0 0 0 v * 

x sin(/z)e-"' ,3T sin(7ny)e-"'' m2T {sin(na;oi) - sin(nx 02 )} e-^^dTdndmdl + 

OO OO oo 

* Mw) -» f (fer) -for) -for)}*— 

(vt) Multiple lines of finite lengths [z Q2lt - 2on], [*02t " x on] and [y 02t - y ou ] passing through (x 0t ,y 0 t) for 
i = 1,2...,L, (yot^oc) for c = L + 1,2,..,M, and (x 0t ,2 0 t) for i = M + 1,2...,N respectively. Where 
(L<M<N). 

_ = 1 

*X> <«> {-(«<*) - [ic„y{^£ + + »»*)} - 

.^ { <^ +l ^ }(s+ „^_^ { ^ + ^ }(s+ Jj + 



+ 1 



Af 00 T / / 
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y«^ yz (m, /, a) + my v ^ xz (n, /, s) + £.^ xy ( n , m , a) ) 

(a + ij.n' + t^m' + ^P) >sin(nar) sin(my) cos(Zs)dndmdl + 



+ 4 



oooooof SS 
™7x ' 




\ Tlx T »?y . »7, J 



/i(l/{^ +!E S a:+!£ ^} i ) 



I ^ »7y ^ T7* / 



dudvdw 



(10.4.17) 



and 



"{-(^)-(^r)*-(^Sf)-(^)}* + 



1 M y. Ut 
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« ( t0) + - ( t0) - (^5) - - (35) } * ♦ 

i 1 j J ( i , \ f Qt (t - hi - «) f -Iffifltl! , ,('+'et) a 1 f 

W^jki 0 / — : — l +e ^ 



U / J 



x (erf f + erf f - erf f V -^) - erf ( L+JSS^ \ du + 

I VZv^y \2^uJ \2y/rwi) V2V^u/J 

00 OO 00 t 

+^ // Jj{ n ^v* KM- T )+™?v^*z (n,/,t-r)-f ^-^ XJ/ (n,m,t~r)|x 

0 0 0 0 ' 

x sin(7zx)e-^ n3r sin(m2/)e"^m 3 r C os(Zz)e- ,| ' |ar dr Aidrndl + 

H t- ■ / / / (p (w, v, w) [e~ ( - e" ( lie" * - e~ 1 x 

W^^W^ £ £ I J I J 

e -1- e 4^it ^ dudvdw 

The spatial average pressure response of the line [z 0 2<> - *oi<>]> t = is given by 
1 



(10.4.18) 



^CtTT 2 (ZQ20 - *01o) y/^hVy 



x Y, q L ( s ) e ~ 5t ° 4 / Z2 { sin ( W2 02o) - sin(u2 01 o)} {sin(uz 0 2t) - cos(uz ou )} x 
11 0 

-* ^} <•♦«*.>}+* ♦ hFF^}]- 

+ 1 

x£> (.) e-^|^M {cos(uxou) _ cos(Kx 02t )}/Tir 0 |^{(i^ + k^i)!J (5 + u277x) | + 



(j>Cf7T^ (ZQ20 - *01$) 

N 00 . *02O 



00 2020 r ✓ r- " N 

x£ ft (.) e-f!^ (cos^on) - cosiuy^f L U{ <*^Z + ^f^}( S + «»*)} 

0 *oio ' 
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p = 



+— * 7771 / ™7»^ (m, f , 5) + mi^ az (n, /, s) + jj ,,, ( ra , m , ^ | 

x sin(nx) sin(mj,) {sin(fz 02O ) - sin(/zo 10 )} dndmdl + 

^tt^^}^}-*W fc s £+ H a: }«- + 'w} + 

. /r J // (* ± ") 2 (2/ + ") 2 l 11 

IV t + ~%~r + /2??2) /J ^ sin (^2o) - sin(/zoi^)} didudvdw (10.4.19) 

and 

= 1 

80C t 7T (2Q20 ~ ^01<>) y/Whl * 

/'^^{ e -^- e - !! ^}{ e - ! ^. e -^j )< 

x (erf ( ^Sh) + erf / x+xon \ _ , / x-*o 2t \ /* + x 02t \ 1 J 
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x(erf(^^^_ erf /'£°!±^^ +e ^ f /*>io+*oi\l 

g oo oo oo t 

^(ateo-zoio) I fill { nr, ^v {m ' * " T ) + m ^** (n, I, * - r) + -L^ („, m> * - T )\ x 
oooo 0ct J 

The spatial average pressure response of the line [x 02O - x 01< >), * = <>, is given by 
p= _ 1 

^W{^^}<— >}]— 

i m ~ 

+ teir a (x 02<> - xoio) JJf ( s ) e ~ sto j^{<*s(ux ou> ) - cos( U x 0 2o)Kcos(«x 0 n) - cos( U x 02t )}x 

I \ 

<f>ct* 2 (x 0 20 - x 01 <>) v/^T X 
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111 l / n7?x ^v* *. •) + m ^xz (n, /, s) + (n, m, s) 1 

J J J n| (s + »7»n 2 + V" 2 +^ 2 ) J 



(*02<> - £01$) 

x cos(iz) sin(my) {cos(nx 0 io) - cos(nx 02 o)} dndmdt + 
~ K °\)f{ + } ( S + }] sin ( nu ) {cos(nioio) - cos(m 0 2o)}dn<iudvrf«; (10.4.21) 



and 



P = 



■K^) --( w) --( * 

<Mw)-»<w) + Kw)-Kw)W 



8(^Ct7r (z 0 2O - zoio) V^Ty 
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* Mw) - ( W -(^) - ( W)} * 
*Mw) + -fe)--te)-^)}- + 

g oooooo i 

+ * 3 (x 02<> - Xoi<> )J J J jn \™l3v» ("». U - r) + m^„ (n, /, t - r) + -Lf („, m, * - r)l x 
oooo v 0 * J 

x {cos(nx 0 i 0 ) - cos(nx 02 o)} e~"* n ' T sin^e-*" 1 '' cos(iz) e -'«' : ' T drdndmd/ + 

M^) + -(^) -(W) -(^)}*-* (UUfl, 

10 ' 5 » > 0^iT>0 10 ' 4 ' 6XCePt iDitial PreSSUre P (!B ' S/,Z ' 0) = P " a CODStant for aU x > 0, 
The solution for the continuous point source is 



uus point source is 

to ( 27r ) * y/VxVyVz /l£^xo^ , {y^yo^ , l £= * 2 P\i 
L I J 



I *?v H» J 11« 



(ix+xoli + .ij-a^li H f (»-»„)» 

/Ix+xo^ + (y-y 0? * ^ (z+zoP ji ( (x _ xa)2 



Vx ^ i) v ^ v* J 



*fr ^ *?* 



I »»» 7, ^ f/, J 

+-1/7/3 n ^»« ( m - <• *) + m »?^xz (n, /, a) + ^ xy (n, m, ,) ) 
* 3 {JJ\ (5 + ^2+^2+^/2) Wnx) sin(my) cos(lz)dndmdl + 

{ / »+t|„m» I 
1-e V - I sin (my) 
T _ , / ■ 

» 7 m(s + »7 v 7n2) dm (10.5.1) 
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and 



t-t 0 



^»vwp / n r "^- e ^j* 

x 4 e + e ***** > dr + 

00 00 oo t 

+^3 J J J J yiVz^yz (m l / l *-T)+m^ x , ( n J,*- r ) + -L^ xy (n,m,*-r)}x 

oooo J 
x sin(nx)e- ,?aenaT sin(my)e- T? « m2r cos(/z)e" ,7 « ^drdndmdl + 



2— OO OO = OO OO 

where ^ (n, !,*) = // (x, z, 5) sin (nx) cos (iz) dxdz, ^ (m, Z, 5) = J / ^ (y, z, 5) sin (my) cos (iz) dydz, 

_ 00 00 
= 0000 _ 0000 

^ v (n,m,5) = / fip xy (x i y 1 s)sm(nx)sm(my)dxdy y 1> xz (n,l,t) = / / V*z (x, z, t) sin (nx) cos (Iz) dxdz, 

00 00 
= 00 00 _ 00 00 

*l>y Z (m,l,t) = f Jil> yz {y i z 1 t)sm(my)cos(lz)dydz i (n, m, *) = / / y> xy (x, y, *) sin (nx) sin (my) dxdy, 

00 0 00 
^ (x, z, 5) = Jip x£ (x, z, Q e-"<ft, ii yz (y, z, 5) = / $ yz (y, z, t) e^dt, and ^ xy (x, y, 5) = fx/, xy (x, y, t) e~ st dt. 

0 0 0 

(i) A line of finite length [zo 2 - *oi] passing through (x 0j yo). 

q(s)e 



00 

-ato 



m 2 + 77 /2) >sin(nx) sin(my) cos(i.z)cfoidmcH + 



+ K 0 

8 7771 ( m > l ' a ) + m Vy^xz ( 




OOOOOOf = 

R 

+ 



(5 + 7] x n 2 + 77 y m 2 



ool 1 -^ V' + - m 1 sin (my) 
+ 2P/ 7] J 



?r 7 m(5 + 7/ v m 2 ) 

and 

t-t 0 



n _ ^(*-*o) / <l{t- to- r) f .ujaiii -<-.+«n>M f -(yv<)> 2 -<»+*o) 2 l 
P - q~t : / < e - e e 4r »v T - e 4 "i> T ^ x 
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oo oo oo t 



w uu UU l 

+ ^ / // / { nvx * vz K 1,1 ' r) + (n - M " T) + (n - m ' * - *■>} x 

x sin(nx)e-«'»°T sin^e""^ cos^Je-^'drdndmd/ + 

+ P /erf(-^erf(-JL- ^ 

\2^<; V2^/ (10.5.4) 

The spatial average pressure response of the line [zfa - h obtained by a ^ .^.^ 
(7 (5) c — y^j r f / "" 



+ K 0 ( 

8 



7 7 7 1 ( "^ ^ «> + (n, /, s) + ^ xy ( n , m, s) } 

{ { { 1 1 l» + «h» i, + vn ,, + ifc/ 3 ) J 



(*02 ~ .201) 

x sin(nx) sin(my) {s in(/^ 2 ) - sin (/z 01 )} dndmdl + 

f - X J 1 

ooW-e V " >sin(mj/) 

+ ^ /i / ■ 

* { m( s + Vym 2) am (1055) 

and 

P = ^r^T^)^/ 5 — ^ [(^2-^)erf (^f^i) - (z 02 + ^erf (^=^) + 

OO OO OO t 

+ t 3 (*02-zoi) //// 7 { n77 *^ (m ' < ~ r ) + (n, /, t - r) + -LZ („, m , t - r) 1 x 

0000 <PQ 'J 

x sm(nx) e -^ si„ We -^ {sin(/2o2 ) _ e -,^ dTrfndm<fl + 

+p,erf( * W-W) 

\2vW V2 V ^ > / (IO.5.6) 

(»). A line of finite length [x 02 - x 01 ) passing through (y Q , *,). 

OO 



82 



W y y i (3+» ? xn2+ J?v m2+^ 2 ) r in(nx) sin(m!/) cos(z ^ dridrmi/ + 

— , 1 , 



* J m(s + Vy mi) dm (10-5.7) 

and 



oo oo oo t 

+ ** 0 o oo { nr)X ^ VZ *' * ~ T) + m7?v ^ 12 (n ' * " r) + l^*v {n ' m ' 1 ~ r) } x 
x sin(nx) e -"* n:,T sin(mj/) e -">' mar cos(Z2) e -"«' : ' T drdndmdZ + 

^(s^MsTp) (10-5.8) 
The spatial average pressure response of the line [z 0 2 - x 0 i] is obtained by a further integration. 

-*{^ + ^_}c + ^}]* + 

, 8 7 7 7 1 J ™7»^,* (rn, f, 5) + m77 y ^„ (n, /, s) + ^ ay (n, m, s) } 

n 3 (x 0 2-x 01 )J J J n| (« + ffcn a +V* a +'?»f a ) ? x 



x cos(te) sin(my) {cos(nx 0 i) - cos(nx 02 )} dndmdl + 



+PI 



i f 1 - e-»Vt) 2 v^ / i_ J _ . 

•V y t^o-xoio)/ m(. +V »a)l *" 



(10.5.9) 
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and 
P = 



4<f>Ctir{ 



w^-^-te?) - •) — m — (^) - 

oo oo oo t 

+ **(xm-zoi) JJJJn \ nrtx ^y* {m > " T > + mT >v$*z (», i, * - r) + ±^ xy (n, m,t-r))x 
0 0 0 0 w J 

x sin(/z)e-"«' JT sin(my)e-"«' m ' r {sin(m 0 i) - sin(nx 02 )} e-* B ' T <*rdiu6iw« + 



+Pj erf 



^020 er f 



(10.5.10) 



(its) Multiple lines of finite lengths [z 0 2* - *on], [*02i - *on] and [y 02t - Von) passing through (x 0t , y 0 .) for 
i - 1,2...,L (y 0 i,^0t) for t = L + 1,2...,M, and (xom^i) for ^ = M + 1,2..., TV respectively. Where 
(L < M < TV). 

1 



V = 



<$>W 2 yJrWh 



fat* 2 y/TlxVy 
N 



' (x-so t ) 2 (z- 
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(s + wi + rtymt+vJ*) >sin(nx) sin(my) cos(lz)dndmdl + 

00 \i-e~ X yr^\ sin (my) 

.2P/ f { J 

W m( 3 + t ,X) ( 10 - 511 ) 



000000/* — — 

& [ f f) n7 lx^ yz (m, /, 5) + mfy^ 



0 

and 

t-t 0 

p = 



8^*vWS 7 « I ^^- e ^^j|c ^^^-e 

1 M *~*o* 

8^^v^ t= 4ti / » I 6 J I 6 " 6 j x 

"{-(^) + -(^)-(^)-(^)}*. + . 

x /erf f ^Mi'l + erf ( y + yo h\ _ erf ( v~V^ \ _ erf f lf + ltoA 1 . . 

^ 00 00 00 t 

J J /{ n ^^( m ' Z - < -- r ) + ^y?xz(n.',<-r) + -i-^ xl/ („,m,i-T)}x 
0000 9°* J 

x sin(nx)e-"' nJT sin(my) e -">' mJT cos^e-"' 1 ** drdndmdl + 

+ ^ CTf (2^)^(2^) (10.5.12) 
The spatial average pressure response of the line \zo 2<> - zo X<> ], t = 0, is given by 
1 

P = x 

4>cfir 2 (zo 2<> - ^oio) y/fh% 

L 00 

x £ * (a) e"" 0 ' / -j {sin(«zo2<>) - sin(uzoi<0} {sin(«zo 2l ) - cos(u^n)} x 



85 



X 



-f \ x 

<t><hn 2 (2020 - 2010) V^Mfo 

+ 1 

>}]— 

7 7 7 1 l 7 "?*^/* K U) + (n, t, s) + ±$ xy (n, m, s) ) 

I I I 1 \ (s + T, x n* + Tj v m*+T, z P) J x 



(2020 - zoio) 

x sin(nx) sin(my) {s\n(lzo 2<> ) - sin(/«oio)} dndmdl + 
1-e V - }sin(my) 



- Q m(s + 7? y m 2 ) am (10.5.13) 



and 

1 



80Ct7T (2020 - ^Oio) y/V^h 
X * 



" " erf (^) - (-0 - *o erf (a^p) + 

+2^|«- ( *Uy>* _ e- 1 ^ 21 ^^ -e-^^ 22 ^ 2 ^ + e- ^^2i ^ ! J _ 
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• KV) ♦-( W) - -(W)} - 

-W^») + -(?5»)-gai)-(^)}* + 

g oo oo oo t 

+ «>(zo 2 o-z 01<> ) J J III {m ' l > t ~ T ) + (n, Z, t - r) + ±^ xy ( n , m , t _ r) | x 

x sinMe— ^„)r^ {s in(/ W - sin(^) } e'^rd^ + 

+P/erf(-4-)erff-X_) 

\2vW V2vW (10.5.14) 

The spatial average pressure response of the line [*„,<> - x 010 ], . = <>, is given by 



P = 1 



i w 7i 

+< ^^ (*»«0 - *dio) VW t= 5?' (5)e " 7 ^ {cos ^ UXo ^) - cos(uio2 0 )}{cos(ux ou ) - cos(«xo2 t )}x 
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^tt 2 (a: 0 20 - xoio) s Jryru 



+^7—? , 7 7 J 1 / "fc^("*.'.«) + "^?..(n,l,«) + (n,m, a )) 

* (X020 -x 01<> ) / / y n j i»+n.f*+v»'+n.z») J x 

x 008(b) sin(my) {cos(nx 010 ) - cos(nx 0 2o)} <*ndmd/ + 

sV y ^(^02o-x 01<> )y TTrm^i dm 



m(s4-%m 2 )^ 

and ( 10 - 515 > 

- WW) - ( W) -'-(W) -OWN - 
"W^M^M^)--(^)}-. + 

_j_ ^ 

80C t 7T (X 0 20 - Xoio) yfifcTfo X 
Af *~*o* 

x „?/ (, - W / Sj£ ^^{'-^ + «-^}{«- tt 8» e -«- ft «* l! lx 
+ «*»+*)-( a ^)-**. + *0-f(^s) + 
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* {- (W) - - ( - - ♦ - (^) } x 



oooooo t 



* 3 (xo 2 o-soio)////n (m> * " T ) + (n, I, * - t) + -L^ B (n, m, t - r)}) 

0 0 0 0 J 



x {cos(nxoio) ~ cos(ni 0 2o)} e~ r, * n ' T sin(roy) e -' I * mar cos^e""' 1 ' 7 drdndmdl + 

< e — e > 



+P/ erf 



(2^) 



*-(55)-*-(5&) 



2s/nr) x t 



(10.5.16) 



10.6 The problem of 10.1, except for all t > 0, p(0,y,z,t) = * v , (y,z t t), y > 0, z > 0, 
p(x,0,z,t) = (x,z,t), x > 0, z > 0, and Msg&H -A xy p(*,y,0, i) = - fx « t) 
x > 0, „ > 0. The initial pressure p (., y, 0) y , .^(x.l, z) ind Its &rttive Jend 
to zero as x — ♦ 00, y — ► 00 and z — » 00. 

10.7 The problem of 10.6, except the initial pressure p (*, y, z, 0) = P/ , a constant for all x > 0, 
y > 0, and z > 0. 
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The invention being thus described, it will be obvious that the same may be varied in 
many ways. Such variations are not to be regarded as a departure from the spirit and 
scope of the invention, and all such modifications as would be obvious to one skilled in 
the art are intended to be included within the scope of the following claims. 
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